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ABSTRACT 


This thesis is an investigation of the earth's 
hatvural electromagnetic field in the frequency range .O01 
CO 5 cps, ibs application to devermine the conductivity of 
the Garth at depth, and a study of spectral techniques as 


a method of analysis. 


Magnetotelluric results from central Alberta give 
crustal resistivities of about 1000 ohm-m with no indication 
of a.major decrease in, resistivity abt, the Mohorovicic dis- 
continuity (or at shallower depths). In southern Alberta 
Pie orupbeal Teal SLiyities are also found to be about 
1000 ohm-m to a depth of 35 km where there is a decrease to 


30 ohm-m. 


Data from a single station near Montreal Lake in 
central Saskatchewan have been interpreted as indicating 
that the station is on the conductive (60 ohm-m) side 
of a vertical fault striking NNE and extending to a depth 
Ore i0eto 25 kn. The resistivity contrast across the fault 
may be of the order of 100/1. The basement is more 


conductive than either surface layer. 


A study of the large vertical magnetic field at 
Montreal Lake indicates that sources are important for 
individual pulsations. However, averaged over long sections 


of record, the geologic effects may predominate. 
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For Montreal Lake records with Ky = 3+ or less, 
the energy falls off quite rapidly up to 0.14 cps and then 
decreases very slowly over the frequency range to 5 cps. 
The spectra with larger Ky decrease much less rapidly to 


BHOUGs Ws eepe OULD tO a -lower energy level. 


Digital spectral analysis was found to be a 
satisfactory quantitative tool. Prewhitening is effective 
in improving estimates but trend removal must be used with 
caution. PManbelbiestimatesirare more satisfactory: than 
hanning as all the estimates are positive and the 


coherencies are more meaningful in the case of noisy records. 


Coherency was found to be only a rough measure of 
the quality of the data as even for very large values 
(> .99), the orthogonal electric-magnetic ratios are not 


necessarily smooth. 
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Spec hrm of . square wave of frequency 
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Illustration of HF spectra without 
prewhitening (CS TUKS) 
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Tllustration of HF spectra with 
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CS TUKS coherency as a measure of 
spectral quality 
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Chapter 1. INTRODUCTION 


This thesis is an investigation of the earth's 
electromagnetic itield.in the frequency range .0O01 to 5 cps 
(the region of micropulsations and their induced telluric 
currents) and its application to determine the resistivity 


ef une earth as a function of depth. 


Micropulsations and Telluric Currents 


Micropulsavions are Tluctuations inthe earth's 
magnetic field in the frequency range .O001 to 10 cps 
with amplitudes varying from tens of gammas (1 gamma = 
107°? oersted) at the low frequencies to .001 gammas at the 
nigin-frequencies. While interest “In this frequency région 
has blossomed over the past ten years, investigations in 
most ranges were first reported over twenty-five years ago. 
Stewart (1861) and Eshenhagen (1897) reported pulsations 
with 10 to 50 second periods; Ebert (1907) and Sucksdorff 
(1936) observed oscillations in the 0.1 and 2 to 3 second 
period ranges respectively; while Rolf (1931) reported 
the giant pulsations in northern latitudes whose periods 


exceed 1 minute. 


From experimental knowledge it is recognized that 
micropulsations can be divided into two classes: those of 
a regular and mainly continuous character and those of an 
irregular or transient character. The most recent 


recommendation to the International Association of Geomagnetism 
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and Aeronomy (Jacobs, Kato, Matsushita, and Troitskaya, 
1964) proposes that they be designated Pe and Pi 
respectively. » Itvis further«proposed that they be 
divided into period ranges based on their physical and 


morphological properties. Tables 1 and 2 show these 


classifications. 
Lawe.. Table 2 
Notation for Regular Notation for Irregular 
Pulsations Pulsations 


Notation Period Range, sec Notation Period Range, sec 


Pe 1 Qo2-5 Pil 1-40 
Pe 2 5-10 Pi 2 40-150 
Pe 3 10-45 
Pe 4 45-150 
Pe 5 150-600 


It should be pointed out that in any such classification 
Gaeenoundary divisions are rather arbitrary. Further, 
pulsations of different types often occur together. For 
example, pulsations during a magnetic bay (Pb) can be 
expressed as 


Pie ete PC. oe 


Although these classifications are convenient for observa- 
tories, the analysis techniques employed in this thesis do 


not lend themselves to these divisions. 
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In one of the earliest theoretical papers 
Stormer (1906) suggested that the charged particles 
contained in a high speed cloud emitted by the sun are 
captured by the earth's main magnetic field and spiral 
about the earth's magnetic line of force. Chapman and 
Ferraro pointed out that the charges would repel each 
other and the streams would disperse. Stormer's 


theory has been revitalized by the work of Bennett and 


Hulbert (1954) on magnetic focusing. Recently workers have 


used variations of this model to explain pulsations of 
type Pe 1. Wentworth and Tepley (1962) have suggested a 
fast electron bunch model while Heacock (1963) suggests 
a fast proton bunch model. In both cases the emission 
frequency is associated with the frequency at which the 
bunch oscillates between magnetically conjugate mirror 
points. The conjugate point phase shift should then be 
close to 0°. Tepley (1964) however reports observations 
with a phase shift of 180°. Jacobs and Watanabe (1963) 
proposed that the emission frequency is associated with 
the resonance frequency of hydromagnetic waves generated 
in the lower exosphere with the resonance oscillations 


excited by proton bunches oscillating slowly between 


magnetically conjugate mirror points. In this case a 180° 


fine-structure phase shift between conjugate points is 
predicted. Unfortunately this model predicts that the 
emission frequency should increase with geomagnetic 


latitude which conflicts with experimental observation. 
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Holmberg (1951) proposed that the sources of pulsa- 
tions are intralayer hydromagnetic oscillations within 
the E and F layers of) the.ionospheres..-.Although. this 
particular trapped mode has been rejected because of the 
high attenuation of hydromagnetic waves in these regions 
(Piddington, 1959), hydromagnetic waves continue to form 


the basis of most theoretical work in this field. 


Dungey (1958) considers hydromagnetic waves in 
the magnetosphere and has concluded that two modes of 
oscillation are possible, (1) toroidal oscillations 
with each surface of revolution of a geomagnetic line of 
force oscillating independently, and (2) poloidal oscilla- 
tions.with the magnetic field oscillating, in. the,meridian 
plane. _He suggests that these disturbances may. be 
excited by the flow of interplanetary. gas, over. the surface 
of- the magnetosphere. Because of,the complexity of the 
equations in spherical coordinates, the coupling terms are 
usually neglected. -In theory, the periods of the 
fundamental and higher modes may then be calculated and 
compared with observations. Unfortunately, the Alfvén 
wave velocity in the exosphere is not yet well enough known 


to do this with great accuracy. 


Watanabe (1961) associates the poloidal 
oscillations with Pc and the toroidal oscillations with 
the giant pulsations. He suggests that for non- 
axisymmetric toroidal oscillations, radial oscillations may 


be induced which are observed as Pi. 
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MacDonald (1961) points out that the sharp 
maximum in the Alfvén wave velocity at an altitude of 
3000 km may allow the trapping of energy between this 
maximum and the FP, layer of the ionosphere. This allows 


Baurcer vertod OBCLILatLlone or "types Pa land *P¢e 2. 


Jacobs and Watanabe (1962) postulate that the 
hydromagnetic waves generated in outer space are propagated 
along the magnetic lines of force towards the polar 
regions where the lower exosphere acts as an amplifier 
at selected frequencies. Secondary waves are then 
emitted from the regions of incidence and propagated to 
lower latitudes... Their calculations show that periods 
corresponding to Pc 1 may be amplified by this machanism. 
Their suggestion that type Pe 3 might also be explained 
by this process has been supported by the more 


detailed calculations of Prince and Bostick (1964). 


The electric field amplitudes range from a few tens 
of mv/km at low frequencies to .001 mv/km at high frequencies. 
In 1921, Davy postulated the existence of the induced 
telluric currents. Barlow (1849), however, first obtained 
convincing evidence that such currents flow in the earth by 
studying the disturbances created in telegraph lines. He found 


that these disturbances are most pronounced during periods 

































qtste eft tedt tuo sdntodg (2307) SienoteaM 

46 spud itis np ts yiloolsy svaw nawILA sts w 
elie neewted yaitene to patqasid stiy wolle ¥ 
awolle gict .sterkeonet arig ‘to se yet ,t elt 
© of bas (09 Beqys to enotteilioas bobreg 





sit dao¢ stsivteoq (Sd0L) sdensten Das adoogt 
bevarsqoiq srs so9sqe ‘= evvoO ait Betataisa sesvew of 
asiow sid ebxswot gore? ‘io eemtl obtedgam sag 
qTSstitfame ns es eton stedqeoxs tewol ais 

nent Sis eevew yasbwdsse .selsreupst 

ot bevsusqota Deis sonsbiout "te eroizet edt 
abotveq: dsdd wore snoltsivoias sient? .at 
smatoadosm etde yo Delt ifems od yew L of oF 
poctslaxs sd oasis Soain € oF egy’ ends not 

stom sfi Yd bedroagiva onded aad 

.(#80L) oldacd base Soniva to enottatac 

a” 

ened wet s mort sgaet ashuttlame bislt oftdosls oe 
eetomeupet? dgin ts mivm 100, of solonanpett Wot. i: 
besvbnt sid to sasatetne odd beveled. ee ad. 
benietdo satit vrevewor .(@H6i) woltes 4 rere) 

yd cifase sft ai wolt sinerivo dove tady e568 
Bauer oH .eentl dgetgelet af pedaers rr 


abolteq ee ete | 
- ie : 
i hale 

ed o> 7 





4 
vis 


“on 7 
a yp ois 








of magnetic activity. By the end of the nineteenth century 
many observatories recorded telluric currents. While the 
correlation of telluric currents and micropulsations was 
appreciated, thesdifficulties of analysis and instrumentation, 
coupled with a seeming lack of positive results, caused a 
Pinictvuation sof interest in the subject .and consequently «in 


thewnunber .of recording ustations : 


The Telluric Method 


The first .application of these ULF electromagnetic 
fields to determine the electrical characteristics, of the 
earth was the telluric method proposed by Schlumberger in 
totes elie Waste Of CHLS bechnidue 1s that the tertluric 
currents will flow around areas of high resistivity. The 
quantity measured is the change with time of the average 
electric field between a pair of grounded electrodes a few 
hundred feet apart. At a base station two measurements at 
right angles are made. Simultaneous measurements are made 
with identical arrangements at various locations in the 
area of interest. Current patterns may then be mapped and 
large conductivity contrasts such as highly resistive salt 


domes located. 
The Magnetotelluric Method 


The magnetotelluric method uses the horizontal 
components of both the electric and magnetic fields at the 


surface of the earth to determine the electrical resistivity 
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as a function of depth. Tikhonov (1950) suggested the method 
as a tool for exploration to great depths while almost 
simultaneously Kato and Kikuchi (1950) showed that measure- 
ments made on phase angles between the electric and magnetic 
fields could be explained on the basis of a two-layered 
earth. In subsequent papers (Rikitake, 1950; Tikhonov and 
Lipskaya, 1952; Lipskaya, 1953) the original theory and 
interpretation was extended. However it was not until 
Cagniard (1953) published his comprehensive paper on the 
theory and graphical interpretation techniques that the 


potential of the magnetotelluric method was realized. 


In the basic magnetotelluric theory as proposed 
by Cagniard, the two fundamental assumptions are that the 
earth is horizontally layered with each layer being 
electrically homogeneous and isotropic, and that micro- 
pulsations are plane electromagnetic waves impinging on 
the earth. Under these conditions Cagniard showed (see 
Appendix A) that the resistivity may be found as a function 
of depth if we know the apparent resistivity 

a 2 
p, = 0.2ETIET/IEI 
Beee UNCC LOM, oO. period, where 
Pp, = apparent resistivity (ohm-m), 


T period (sec), 


lI 


II 


horizontal electric-field amplitude (mv/km) at 
period T;, 


| H | = horizontal magnetic-field amplitude (y) at 
period T. 
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He also showed that the phase difference between 
E and Asi eouid- pe usped*for-interpretation. In this paper 
Cagniard developed the theory for l1- and 2-layer cases. 
Berdichevskii (1960) extended the formulas to n-layers 
while Yungul (1961) has presented 3-layer interpretation 
curves. Solucions for the effects of a fault on the earth's 
electromagnetic field have been presented by Neves (1957) 
and d'Erceville and Kunetz (1962). Rankin (1962) has 
extended the theory to the case of a dike. These results 
may also be presented in terms of apparent resistivity. 
For application of the magnetotelluric theory it is then 
necessary to Tilnd E and H asses .LuUncliam of period, Irom 
which the observed apparent resistivity may be calculated 
and compared with theoretical models. -Neves points out 
that for two dimensional structures additional information 


is contained in the phase. 


The fundamental advantage of the telluric and 
magnetotelluric methods over conventional electrical pros- 
pecting methods is their depth of penetration. The conven- 
tional techniques are limited by the inverse square law 
for point sources, whereas the penetration of plane 
electromagnetic waves in a uniform conductor is given by 
the skin depth 


1 
ea (Fa.0 owl 
km). 


where p = penetration ( For ULF electromagnetic 


waves and typical resistivities (1 ohm-m to 10° ohm-m), the 
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penetration is many kilometers. A second advantage is that 


nature provides the signal. 


Wait (1954) first investigated the effect of 
non-plane wave sources on the magnetotelluric theory and 
found that the correction terms contained second and higher 
order Spatial derivatives. In the case of oscillations of 
10 sec period over a homogeneous earth of resistivity 
1000 ohm-m, he found that these terms become important if 
E and H vary appreciably in a distance of 35 km. He 
pointed out that if the micropulsation sources were in the 
lower ionosphere at an altitude of 100 km, corrections 
would be important for periods greater than 10 sec. On 
the basis of observational evidence, Cagniard (1954) 
insisted that the fields are uniform over hundreds of 


kilometers. 


Price (1962) investigated the effect of the source 
dimension on the intrinsic impedance | E,/Hy ly Mor @ 
homogeneous earth of conductivity o, he finds that 


(MKS units) 
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where v= separation constant. 
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10. 


He points out that 27/v can be taken as the linear 
dimensions of the source field. For v = 0, the source is 
infinitely large corresponding to Cagniard's plane wave 
solution. Wait's result is obtained by taking ve<cws. 

To obtain an estimate of the smallest and largest values of 
v which can occur he equates 27/v to the circumference 

of the earth and to 4 times the height of the ionospheric 
currents. He concludes that the values of v of interest 
will lie in the range 1.57 x 107 km to 1.57 x 107 um. 
The calculations of Quon (1963) using dipole sources confirms 
this estimate for the upper limit of V. Price then 
considers a two-layer earth with a conductive upper layer 
and Intinively resistive substratum. In this situation he 
shows that Wait's estimate of the source effect is too 


small and that even for sources of global dimensions 


correction terms may be required at long periods. 


Since. the publi cathionsot. CagniLaerd)s classic 
paper, the results of more than 25 magnetotelluric 
soundings have appeared in the literature. Because of 
adifficulties in detecting the low amplitude signal at higher 
frequencies, most investigators have confined themselves 
to the frequencies below O.1 cps which correspond to depths 
where relatively little is known about the resistivity 
of the earth. Due to the possibility of source effects 
at these periods and the wide scatter of values obtained in 
a single experiment, there remains some uncertainty as to 


the validity of the interpretations. 
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In addition to the difficulties of signal detection, 
the problem of quantitatively analyzing the data has been 
considerable. In almost all cases, workers in this field 
Have UBea- Charu recorucrs. “lr order to obtain sufficient 
accuracy at all frequencies, this requires that the response 
of the detecting system must be closely matched to the 
character of the signal; otherwise, some part of the spectrum 
is lost because it is either too large or too small to be 


read accurately from the chart. 


While the technique of picking sinusoidal events 
from records has been widely used (Scholte and Veldkamp, 1955; 
Niblett and Sayn-Wittgenstein, 1960; Wiese, 1962), it has not 
proved very satisfactory. For one thing, the period is 
never accurately measurable as there is no single period 
to a discrete event. Secondly, .large sections of record 
must be scanned before observing sinusoids over the complete 
frequency range. Finally, apparent resistivities calculated 
using large sinusoids do not yield a smooth curve but 
exhibit a scatter whose range may be up to a factor of 3. 
From this, it is questionable whether the Cagniard theory 


should be applied to single events. 


Other groups have digitized the records and 
found ~ and Hq as a function of frequency using mathematical 


techniques. Garland and Webster (1960) used Fourier analysis. 
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Cantwell (1960) considered the magnetotelluric field 
components as samples of stationary stochastic processes 
and applied the autocorrelation methods of Blackman and 
Tukey (1958) supplemented by cross-correlation techniques 
(Robinson, 1954) in order that coherencies between electric 
and magnetic signals could be found. He assumes low 
coherencies are caused by high instrumental noise levels. 
Using only records with high coherence, it was found that 
comparatively smooth resistivity curves were obtained 

and estimates from different records were relatively 
consistent. He also outlines procedures for determining 
the effects of anisotropy and two-dimensional 
inhomogeneities in the earth. The Electrical Engineering 
Research Laboratory at the University of Texas follows 
similar statistical procedures using analog techniques 
(Smith, Provazek, and Bostick, 1961). Their use of 
magnetic tape recording has increased the dynamic range and 
the elimination of the hand-digitizing procedure has 


removed one source of noise. 


The Vertical Magnetic Field 


In Appendix A it has been shown for plane waves 
incident on a homogeneous half-space that the vertical 
magnetic field is negligible compared with the horizontal 
fields and decrease with increasing period. However experi- 
mental results show that Hy is sometimes as large as the 
horizontal field. For this reason workers have investigated 


H, from several viewpoints. 
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Scholte and Veldkamp (1955) theoretically 
considered reflection of a general electromagnetic field by 
a sphere. They found Ho to be much larger than predicted by 
Che wlayered half-space model but stilw smaller itham ithe 
observed results. In this case Ho increases with increasing 
period and the phase shift between E and Ho depends on the 
form of the field. (For a plane wave impinging ona 


layered earth the phase difference is 0°.) 


The failure to explain the magnitude of i using 
homogeneous or layered earth models has led to investigation 
of two-dimensional structures and particular source effects. Because 
of the sharp resistivity contrast, the sea-land interface 
has attracted the attention of some workers. RokityanskiY 
(1963) has measured the vertical and horizontal magnetic 
fields on the landward side at Mirny. He finds that | Z/H | 
decreases with increasing period for T > 100 sec. Using 
experimental results and the theoretical results of Berdicheskif 
(1961) and Ponomarev (1960) for the effect of wedges on 
electromagnetic fields, he has made qualitative plots of 
the field components near the sea-land boundary. He suggests 
that He should be quite large on the land side but negligible 
on the ocean side. Weaver (1963) attempts to solve the 
same problem using a fault of infinite extent. He finds 
that | Z/H| falls off more rapidly on the seaward side but 
may be as large as the horizontal fields close to the 
boundary. His calculations at two periods, T = 1 and 10 sec, 


show | Z/H | to be much larger for the longer period. 
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Law and Fannin (1961) have considered a line 
source at lonospheric height over an earth of resistivity 
8000 ohm-meters and with a period of 21 sec. The fields 
Nave Deen calculated at the surface of the earth for 
VErl CUSie@ eranrestaromelLhesource .-8 Lt asefeund that He 
may be as large as 50% of the horizontal field and that 
the phase’ between He and the electric field remains 


close to oat 


Quon (1963) has calculated the fields over a two- 
layered Garth resulting from dipoles at ionospheric 
hneizhncs. He finds that He is comparable with the horizontal 
16s a5 the longer periods. The phase relationships are 
Teuner complicaped depending on the type of dipole, 


G@ievance, and period. 


Ob Jjecvives of the University of Alberta Group 


The experimental objectives of the University 


of Alberta magnetotelluric group have been as follows: 


(1) to develop a mobile recording system capable of 
devecting micropulsations and telluric currents in the 


Trequency range .0O01 to 10 cps. 
(2) to develop low noise level analysis techniques, 


(3) to use the system to test Cagniard's theory and 


study the sources. 
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Purpose of this Thesis 


This investigation is limited to the following 


problems: 


(1) the development of analysis techniques 

(2) a-study of source effects as observed at a single 
Station’ on the’éarth's surface 

(3) an interpretation of the field data using the 


magnetotelluric method. 
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16. 
Chapter II. 


DETECTION, ANALYSIS, AND PRELIMINARY RESULTS 


Signal Detection 


Two perpendicular components each of the electric 
and magnetic fields are measured at each location. In 


addition, the vertical H component is recorded. 


Figure 1 is a schematic representation of one 
electric and one magnetic detection system. To find E, 
the voltage difference is measured between two copper electrodes 
approximately 1,000 feet apart. This voltage is fed through 
a i1xediband pass filter to a dc post-amplifier with gains 
up to 30 db and a high frequency cutoff at 8 cps, and 
then to the seven-channel FM tape recorder. Normal 
amplifier settings have been in the range +40 db to +54 db 
in order to achieve the signal level required by the tape 


recorder ( + 1.5 volts). 


The magnetic detectors are induction coils 
eonsisting—of—15,000 turns of wire wound-on.a-permeable 
rod core. The coils are 5 feet long and 4 inches in 
diameter. The submicrovolt signal which emerges from the 
coil is amplified by a split-beam photocell, feedback 
amplifier. This signal is further amplified and low-pass 


filtered before tape recording. 
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Relative time pulses from a chronometer are 
recorded on a direct record channel. To obtain absolute 
time, these pulses are compared with WWV on a high-speed 
paper monitor when reception permits. This allows 


comparison of our signal with those at other stations. 


During recording, two channels are monitored on a 
chart recorder. Using this, the amplifier settings may be 
changed in order that the full dynamic range of the tape is 
utilized and adjustments made for amplifier drift. If the 
Signal level is too low or the bandwidth too narrow to be 
useful in calculating an apparent resistivity curve, the 


tape may then be erased. 


Calibration of the entire system at about 10 
frequencies in the range of interest is recorded at least 
once on each tape (4-12 hours). In the earlier systems 
this was required because of instability in gain. However, 
Sains an the 1962 and 1963 systems were found to be stable 
within a few percent over an interval of several months. 

At each calibration frequency, the output of an 

oscillator (Figure 2) is recorded. At the same time an 
attenuated output is fed through the e electric systems and 
through solenoids, wound around the detection coils, which 
are in series with a precision resistor. From the playback 
of these recordings, phase shifts and amplitude calibration 
curves in terme of Va ae for the electric systems and 


Vout Min for the magnetic systems are calculated. 


Lo. 


cs : 7 
































OL 
Sis Tesemonosto s moxt eeefug enti maces , 
sJifosda aisido of .féringdo badser goetth 8 fo 6 
pesge-dgirl s no VW dtiw betsqmos Sas essing sai 
awolts athT .etimreq noltsqsoot neiw toxt.te 
svoté¢sda tordo ts s8ond riviw Isdgie tve To fom, 
& no betottinom sts elsnasnio ow’ gn ibroost sntuwd ~~ re 


ag Yam agniides sotitiome sit ,eatitd anteU .sebtose sab 
at sgad oad¢ to egcet olmsnyh ([iut sdv seds tebae afb 
gait oT .dttieb solrttiqnuse sot sham atasmgantbs bas be ia 
94 of worrsna ood débiwbasd sd? 10 wel eed af tsvel. 
edd .evtus viivitetsst Instsqqs As sabcesuntan ail 
.beseats od nedd wi 


Of sede ti meteye ottsas sdz to COR oe 
tesel 3s abbxene el sasis3ar to sgnset end af as 
amodgeve tetitss og al .(eavor Sf-4) eqed dase 
198v 90H -itsy at yttitdsetant to sessed bertuper ¢ 
s{dste ed of bouc't stew amedaye FOC hos Ses ortg me 
.adidonom Isisvee ‘to Lsyiewnt as tevo jnooreg wat ot 


as lo duadue sdt .vonaupett no tnemeltiia doe 


oo 


os omit emse edz 3A .bsbroos x ei (S 9% “> i Pes 


an 
a 


_— 


bus emsdeve oitdosis S ant davondd bet at ratio 


cho triw Ailoso noltosdeb sdd bavets- Bovew «eh ii. tc a | 
7 a ; _ 


dosdysiq sit moti .doveteet notetoetq 6. nes 

notisudiiss ebusiiaqns bons attire sesitq . sad ae ” 
bas amedaya atitoels sqt rot. or Ns <j 

-beisftwolso sta amsgeye ‘otis 





Oscillator 






Attenuator 


Attenuator 


TI 


Q) 


Solenoid 






Band Pass 














wound around 





n SHH \ 
SST CEs 


S 


(os 


a 
te! 












Filter induction coil 





s 








\) 


DC Galvanometer 


Amplifier Amplifier 





DC 






DC 





Magnetic 





Amplifier Amplifier 





rape 







Recorder 


Figure 2. Block diagram of field calibration system. 
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For an absolute calibration of a magnetic system, 


the solenoid wound on the coil must be calibrated. This 
is done by placing the coil at the centre of a solenoid 
10 ft. in length, one foot in diameter, and with one turn 
per centimeter. By comparing the outputs of the magnetic 
system when currents are passed through the 2 solenoids 
and calculating the magnetic induction in the large 
solenoid, the quantity Vian ea for the solenoid wound 


SiecrerilieuctLon coll is found. <A calibration curve of 


naa em may then be calculated. 


The 1962 amplitude and phase calibration curves 


are shown in Figure 3 and 4, respectively. 


Basic Analysis Techniques 


The fundamental analysis problem is to find 
pi and H as a function of period (or equivalently 
frequency). Cantwell (1960). and Srivastava et al (1963) 
Gonsidered the,.recorded data as statistical samples, from a 
stationary random process of infinite duration and 
applied the "Tukey" spectral techniques (Appendix B) given 
by Robinson (1954) and Blackman and Tukey (1958). Under 
the assumption that the stationary random process is 
Gaussian and the spectrum reasonably flat, the stability 
associated with each estimate can be expressed in terms of 
the spread in db of an interval containing, with 
prescribed probability, the ratio of true smoothed power 


to estimated smoothed power. For example, the 80 percent 
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range in db is given by 


(80% Lanes ) bee 
af Gb record 1 
(no. of records )[ (ai nation) (resolution) -3] - 


From this the duration of record required for a given 


resolution and accuracy may be calculated. 


It is preferable to consider the data as a unique 
record generated by an aperiodic phenomenon with a continuous 
spectrum. However, analysis techniques for this case have 
not been developed and it must be treated using statistical 
techniques or one of the classical methods where the 
phenomenon is thought of as periodic. In this thesis, 
digital Tukey spectral techniques are used. However under 
these circumstances, the error analysis procedures will 
have little meaning and questions arise concerning the 
required size of the ratio of the number of points In the 
time series (n) to the number of spectral estimates (m). 

In almost all cases presented in this thesis the 


ratio n/m 2 0; 


Since the spectral analysis is carried out on a 
digital computer, the analog data must first be converted 


to digital form. The following analysis procedures are 


followed: 
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24, 
The entire magnetic tape is played back onto a 
6-channel recorder and long sections of record which 
contain wide band signal and are free of transients 
erence ected tor conversion. 
These sections are converted to digital form by playing 
the signal from the magnetic tape through a low pass 
filter (to minimize aliasing problems in spectral 
calculations) and an analog-to-digital converter 
driving a high-speed paper-tape punch (Figure 5). 
During this stage the analog signal is monitored as a 
check on the analog system and conversion operation. 
To facilitate handling, the paper tape is next converted 
to cards using at IBM 1620 computer. During this step 
the paper itape is checked fon parity, Tlagging, record 
length, and magnitude of successive differences with 
the invalid numbers being replaced by interpolated 
values. If the number of errors exceeds approximately 
0.2%, the analog to digital conversion is repeated. 
The data are further edited using the monitor from (2) 
anaa Ulertal plot if necessary. 
The mean and linear trends are then removed and the 
spectra calculated. The spectra, along with the 
calibration data, amplifier settings, and electrode 
spacings are then used to calculate the apparent 
resistivity as a function of period. 
The apparent resistivity curves are interpreted using 


theoretical curves. 
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Preliminary Results~ 


During the 1961 field season, measurements were 
made at four field stations in the Edmonton area - Cooking 
Lake, Onoway, Bon Accord, and Kavanagh (Figure 6). This 
region is a nearly ideal test area for the magnetotelluric 
method. As shown in Figure 7, it is underlain by a great 
thickness of flat-lying conductive sediments, beneath which 
is a resistive basement. There is a vast amount of 
resistivity log information available, including some from 
the basement, as well as gravity and aeromagnetic maps of 
parts of the area. In addition, several crustal thickness 
determinations have been made by seismic methods within 
200 miles of Edmonton (Richards and Walker, 1959; 


Cumming, Garland, and Vozoff, 1962). 


Sinee the conductivity of the sediments. is 
pelieved to bé independent of horizontal direction, it is 
expected that. the. short-period asymptotes of the pair of 
apparent resistivity curves at one station ‘should be the 
same.,, At.the.longer periods the curves may ‘diverge since 
aeromagnetics and basement cores show the Precambrian to 


be grained and variable in composition over short distances 


: The magnetotelluric data used in this section were 
originally processed by the author and interpreted by 
Hasegawa (1962) using regional values for the conductivity 
of the sediments. After a more detailed well-log analysis, 
the data were reinterpreted by the author and published as 
& part of a joint paper by Vozoff, Hasegawa, and Ellis (1963). 
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Figure 6. Locations of 1961 Alberta field stations 
with lithology of Precambrian basement 
as inferred from well samples and gravity 
anomalies (after Garland and Burwash, 1959). 
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Figure 7. Geological cross-section AB (Figure 6) of Edmonton 
area (after Alberta Society of Petroleum Geologists 
regional cross-section GHJ, 1959), with typical 
well-log resistivities for this region. 
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(Garland and Burwash, 1959). In addition to this apparent 
basement anisotropy (due either to anisotropy or to 
inhomogeneity), a thin, very highly conductive basal sandstone, 
whose thickness may reflect the basement topography, is 
present in the eastern portion of the test area. Consequently, 
there are two regions whose properties might cause measured 


apparent resistivity to vary with azimuth at longer periods. 


Fach data channel was digitized separately and 
autospectra calculated with a program using hanning 
windows. The computational definitions used were those 
for the autospectrum in the program CS TUKS (Share No. Blt Vee 
For interpretation purposes, sedimentary resistivities 
were calculated using logs from nearby wells. Using these 
values for the sediments, we then attempted to fit the 
observed curves with theoretical curves by choosing 


appropriate values for basement resistivities. 


Cooking Lake. - The results at Cooking Lake for the two 
orthogonal measurements are shown in Figure ve rhe 
observed curves agree at the shorter periods but diverge 

aa tte Peneer periods. This effect: dan be put into phe 
model by assuming an anisotropic layer exists at the bottom 
of the sediments (2.5 km). Basal sandstone was encountered 
in the nearest basement hole, about 10 miles away, where it 
was 30 m thick and had a resistivity calculated to be 3.5 ohm-m. 


eee 


1 Munk, Snodgrass, and Tucker (1959) outlined these 
computational procedures. 
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This alone is not sufficient to explain the observed effects. 
We must assume either (a) that near the site its actual 
resistivity is lower than 3.5 ohm-m (one ohm-m and less is 

more typical) and its thickness is greater than 30 m in troughs 
elongated in the E-W direction, or (b) that the observed 
aeromagnetic "grain" is accompanied by a similar characteristic 
in basement resistivity. In fact, both factors probably 


contribute. 


Onoway. - For this station it is not possible to obtain 
agreement between the observed and theoretical models for 
either the EL Hy or al HY curve using the measured sedi- 
mentary resistavities for the 2.7 km of sediments and a 


one-layer basement (Figure 9). 


In a further analysis a good fit was obtained 
for the observed he Hi curve by inserting a conductive 
layer at thewbase of the sediments. Since it is very 
doubtful that the conductive basal sandstone extends this 
Targwest,. a conductive upper layer in the Precambrian is 
receired. To 1it the Ee 3 Ay, curve the 12 and 5 ohm-m 
layers must be replaced by a section of 15 ohm-m. An 
extremely resistive basement is also required, assigning a 
high degree of apparent anisotropy to the basement. Such 
anisotropy would undoubtedly affect the shallower apparent 


resistivities. 
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Bon Accord. - As no magnetotelluric information was 33. 
available in the 2 to 10 sec period range (Figure 10), 

an analysis using detailed well-log information was not 
carried out. In fitting models to the observed cair of 
curves it was necessary to use different values for 
resistivity in the top layer. Aside from this complication, 
which is discussed below, good fits were obtained using 

the known basement depth and values for the other parameters 


comparable to those at Cooking Lake and Onoway. 


Kavanagh. - The observed curves at this station (Figure 11) 


cannot be fitted using reasonable parameters. 


Discussion 


It should be emphasized that the Cagniard theory 
does not yield correct interpretations, if resistivities 
depart from the assumed conditions of homogeneous, isotropic, 
horizontal layers. \ However it should still be useful as 


a semi-quantitative guide. 


Interpreted deep resistivities are fairly consistently 
of the order of 1000 ohm-m -- a value which agrees very well 
with those of Niblett and Sayn-Wittgenstein (1960) and 
Srivastava, Douglass, and Ward (1963). By calculating a 
number of theoretical curves with different basement resistivi- 
ties .and.considering the spread of the data, the error in 
interpreted resistivity is estimated to be +250 ohm-m. 
Calculations show that a decrease in resistivity of one 
order of magnitude (to 100 ohm-m) at the Mohorovitié 
discontinuity (40-50 km) is not consistent with the 


observed apparent resistivity curves. However the curves 
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are insensitive to a resistivity increase at these depths. 
Farlier results from the southwestern United States (Vozoff, 
Hasegawa, and Ellis, 1963) indicate an apparent maximum resis- 
tivity of about 300 ohm-m at a few tens of kilometers and a 
rapid decrease beyond that depth. It would seem that the earth 
at these depths is quite different in central Alberta than in 
the Basin and Range area. (Seismic velocities in these regions 
are also quite different (Cumming, Garland, and Vozoff, 1962; 
Eaton, 1963)). The low resistivities found by Wiese (1962) 
from European stations and his summary of the results of others 
confirm our findings that resistivities at deep crustal and 


upper mantle depths vary considerably from location to location. 


As pointed out in Chapter I, Price (1962) showed that 
for a generalized source apparent resistivity curves can depart 
markedly from Cagniard's prediction at "long" periods if the 
resistivity increases with depth. A set of Price curves with 
the geological parameters approximating those of central 
Alberta is shown in Figure 12. For the range of periods in the 
experimental data, it is seen that source effects could only 


: which corresponds to very small 


be detected if v~l07°km™ 
sources. The situation is further complicated as this effect 
could also be interpreted as a conductive lower layer. On only 
two of the curves, Ey -Hy at Cooking Lake and Onoway, is there any 
evidence of such an apparent decrease in resistivity at the long- 
er periods. Since this evidence is based on a total of only 3 


spectral points, and is not supported by the other apparent res- 


istivity curves, it is not thought to be Significant. 
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In a general way these results seem to verify 
the Cagniard theory, at least in the range of periods and 
resistivities involved. However it is clear that the results 
cannot be completely explained by this theory. At the three 
stations where interpretation is possible, the Be as curves 
are essentially the same. However, the BE 7 Hy and 
By.” Hy interpretations do not agree. At Cooking Lake this 
apparent anisotropy can be explained on a geologic basis. 
however at Onoway and Bon Accord, more difficulty is 
encountered as anisotropy must be introduced into the 
sediments. D. Rankin (personal communication) finds similar 
displacements in apparent resistivity curves obtained from 
model experiments. He attributes this to boundary effects 
which would correspond to geologic inhomogeneities. Although 
mile may partlally explain the kesults, it is doubtful 


whether large enough contrasts exist in the sediments to 


explain our results. 


Figure 13 shows three BH apparent resistivity 
curves which were calculated using different sections of 
record from Bon Accord. It is impossible to obtain a 
theoretical Cagniard curve to match the curve obtained 
from Tape No. 3. Furthermore the differences are so large 
that it is doubtful whether they can be explained on a 


geologic basis. 


The discrepancies in the results at Onoway and 
Bon Accord and the rough observed curve at Kavanagh are 


probably due to a combination of factors. 
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(1) At frequencies above 0.1 cps, the noise level was 

often high. In some cases this is probably true at lower 
frequencies. Significantly, our best result (Cooking Lake) 
comes from the record with the largest signal. For very 

long periods some results were unreliable due to drift in 

the magnetic system. It is difficult to say where these 
problems set in. Where the curves became rough or physically 


unreasonable the data was rejected. 


(2) Using hanning, as with almost all spectral windows, 
energy estimates on steep slopes are not accurate. Almost 
invariably magnetotelluric spectra have one or more steep 
Slopes. Again it is difficult to say where on the curve 


thie uUnreliabillity sets in. 


(3) Finite-sized sources may also influence the results. 

Quon (1963) shows that under some conditions the apparent 
resistivity curves which would result if the sources were 
Gipoles might bear no resemblance to Cagniard curves. For 
example, vertical components of H can be appreciaple. 

fie teal LLy 2b 1s expected that the sources of pulsations 
observed at a fixed point on the earth's surface will 

change with time. Thus, the time-averaging power spectrum 
approach is probably averaging fields from a number of sources 


at each frequency. 


(4) At this time there was also some doubt as to the stability 
of the detection system. This was checked by frequent cali- 


bration at discrete frequencies. However, there is some 
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difficulty in picking peak-to-peak amplitudes, since, in 
the magnetic systems, natural signal is always superposed 


on tne calibration signal. 
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Chapter III 


ANALYSIS OF 1962 DATA 


Instrumentation 


From the results of the 1961 experiments it 
was clear that improvements were required in the detection 
system to extend the observable range to both shorter 
and longer periods. The photocell amplifiers were rebuilt 
with an improved feedback system which reduced the drift 
problem at long periods. The use of better dc amplifiers 
Powered cie noige level in the electrics which extended the 
range of usable frequencies as high as 5 eps* «Further, 
the Gainseofsbothrtherehéectriceand magnetic’ systems were 
found to be stable. Thus calibrations were not required as 
often and periods of low signal level could be selected 
for this purpose. Amplitude calibration curves were obtained 
by picking the amplitude of the sinusoids from playbacks 
of the calibration sections in both analog and digital 
form except at very long periods where spectral analysis 
was used because of the poor wave form and low signal-to- 
noise ratio. Phase calibrations have been obtained by 
playing back the calibration sections through a phase 


meter. 


Step off becvbet So liW mMeveve yopabest & 


Sig. 


ans henistk= cobtw gobigosie Sag at Level 


* 


(Fi 
jay stew areltifana Ifssatong od? aber 
















Iti edgar 
agac Sher GO STEWAMA 


eile: ; 


oe E 


” 
oniTequs Loe l env to svivuest ant mort 7 


ab o9fd cl DOtTstesy Galen arcomevotgmet ats 


jrode nied of phat sissvieeds: ome oved 


i voroan 10 2B! act eho sq 


voiced .eqe @ ap dyad se eehireugest 


spo" atie7s76 of Senauin ore oiasosls ote oe a 


ep bax 
i. 
bentley 
BA 
, 





Joatyetha motl ebicaumte ents 45 shuvtiqgms % es 













ttupan Jon wrew ereitasdties asd? ee, 


ch 


atsefee an bingo Isvad teagie wol To abo bt 


avow says obtetdites shyt lig weee 


fedileath bas aelisne cdded of -enaissess 


éviacs let1s9Qs Staaw shoiwet gret 4a 


-o7-fensia vol pep wel gyaw so0q: Sea 


“i benteica Hes even eno idpiitias 


4aafig 5 dauotid ano ijose nots: 


43. 


Purpose and Location of Field Kxperiments 


With an improved detection system it was 
desirable to record data for 4 purposes. 
(1) Data recorded at a fixed location at different times 
were required to test for consistency of apparent 


resistivity estimates. 


(2) Data from a location with good resistivity control 
were required to confirm that experimental results 
could be obtained which are consistent with known 


geology. 


(3) Data were required to investigate the relationship of 


He to the horizontal magnetic and electric fields. 


(4) Data recorded at a site where the sediments are either 
very resistive or very thin, would yleld resistivities 
at upper mantle depths provided that crustal resistivities 


were comparable to those of central Alberta. 


As good resistivity control is obtained by 
using well-logs in thick, conductive sediments, this 
requirement is incompatible with penetration to upper mantle 
depths. For this reason, magnetotelluric measurements were 
made in two areas - at Montreal Lake in central Saskatchewan 
and at Vulcan, Brooks, and Beiseker in southern Alberta 


(Figure 14). 
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At Montreal Lake the well-logs show that the 
sediments are thin and resistive (Figure 15). Since the 
effective resistivity (reciprocal of mean conductivity) 
is 73 ohm-m, frequencies of almost 100 cps would be 
required before this thin sedimentary layer would have a 
marked effect on the magnetotelluric curves. At the 
nearby basement hole, Great Plains Montreal Lake #2-5, 
the well-log indicates a basement resistivity of 200 ohm-m, 
while 13 miles to the northwest at Great Plains McDermott 
West Montreal Lake #5-22, the resistivity is 80 ohm-m. 

The basement rock at these locations is mica schist and 
altered biotite gneiss respectively (Peterman, 1962). 

No other pertinent geophysical or geological information is 
available in the immediate area. However the Canadian 
Shield is exposed 60 miles to the northeast . There, the 
surface geology shows major trends striking NNE. The 
aeromagnetics show similar trends which to some extent is a 
reflection of the presence of pyrrhotite dikes (M. Reford, 
personal communication). The strike of the gravity trends 
to the south of this location is between that of the 


aeromagnetics and true north (Peterman, 1962). 


In southern Alberta, the sediments are thick, 
flat-lying and conductive (Table 3) giving resistivity 


control at the higher frequencies. 
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Figure 15. Sedimentary resistivities at Montreal Lake 
based on well-logs from Great Plains 
Montreal Lake #2-5. 
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Table 3 


Station Information 
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Latitude 
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At all stations the detection systems were 
aligned along geomagnetic north-south and east-west directions. 
Five tapes of data were recorded at the Montreal Lake station 
during the period August 22-24 from which 12 - hours were 
selected for analysis (Table 4). The southern Alberta 
Stations were occupied for shorter periods in early September 
Wea tOLva. of ~ hours being digitized. .In this thesis the 
Montreal Lake data are analyzed in some detail with the 


southern Alberta data (from a-simpler geologic situation 


and lower geomagnetic latitude) being used only for comparison. 


Analysis Procedures 
Data channels were digitized in pairs with 


» H_-E 


4 conversions in most cases -- ok ee 


Hues nC Agen. a. rox 
oe ee youl 
computational economy each section selected for analysis 
was digitized at two rates. The entire section was converted 
at a slow speed for low frequency analysis with subsections 
being converted at a higher speed for high frequency 


Gnaiysis. During this step, the analog signal was filtered 


to minimize aliasing problems. 


After editing, the autospectra, coherence, and 
phase angle were calculated at f = 0, Gs, sell @ alae a 


where — “ he tw the Nyquist frequency. CROSPe (Simpson 


a 
et al , 1961; Galbraith, 1963) which uses Daniell windows, 
was used to calculate all spectra with a number also being 
calculated with CS TUKS for comparison purposes. The 


conversion and analysis parameters for the two types of 


conversion are given in Table 5. 
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Table 4 


Tape Information 


Tape Date, 1962 Recording 3-hr Range Remarks 
Feriod . Jndices ih 
(GMT ) Pp 
Montreal Lake 1 Aug 22 2224-2354 5- 
Montreal Lake 2-1 Aug 23 0556-0727 24, 30 
Montreal Lake 2-2 Aug 23 0751-0843 30 
Montreal Lake 3 (7.55) Aug 23 1355-1554 3+ 
Montreal Lake 3 Marker Aug 23 1608-1747 3+ 
Montreal Lake 4-2 Aug 24 0704-0858 -- Visible 
aurora 
Montreal Lake 4-4 Aug 24 0929-1132 3+ Visible 
aurora 
Montreal Lake 5 Aug 24 1225-1336 dp 
Beiseker 1-2 Sepe 10 0309-0426 Lo 
Velean 2-4 Sent 12 0947-1146 50 
Brooks 1 Sept 13 1130-1343 3=. 55° OT 
Brooks 2-2 Sept 13 1621-1753 U4 


To obtain local time subtract 6 hours at Montreal Lake 
and 7 hours at Alberta stations. 


KX 
The geomagnetic planetary 3-hr-range indices — are after 


Lincoln (1962, 1963). 
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Table 5. Analysis Parameters 


At (spec) n m f,, (cps ) Af (cps) Euany (cps) 
Low frequency 2 1000 100 ek 0.001 UO, 27o 
High frequency 0.050 4000 500 10 0,02 5 


An IBM 836 Document System was then used to plot 
Spectra, coherencies, phase, apparent resistivities and 
spectral ratios. In this chapter these plots (and others) 
are used as examples with a larger collection of the machine 
plotsiappearing in Appendix D. An explanation of the notation 
and conventions» used in plotting is also contained in this 


appendix. 


These plots have been used as a basis for inter- 
pretation of resistivities and vertical magnetic fields in 
Pins Chapter and for a discussion of the effects of 
Simultaneous conversion, trend removal, prewhitening, 


and spectral windows in Appendix C. 


Decimation of Apparent Resistivity Values using Coherency 


Under the assumptions of Cagniard's theory, the 
orthogonal components of E and H _ should be coherent. 
It has been found (Appendix C) that smooth spectral ratios 
and therefore smooth apparent resistivities may be found by 
discarding all points for which coherency-squared is less 


than 0.9 and the points adjacent to these. In practice it 
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is found that this is so restrictive that few points are left 

on some apparent resistivity curves. Since it is preferable 

to have some estimate of resistivity rather than none at 

all, points have been used for which coherency-squared is 
greater than 0.75 or 0.8. This appears to be satisfactory 

in most cases. Before examining the apparent resistivity 
curves, the reasons for variation of coherency and some physical 


implications will be examined. 


In practice, "noise" may be introduced by several 
sources.|.tn any digital record section, noise due to detection 
systems, recording, and analog-to-digital conversion is 
always present at some level. At frequencies where the 
Signal predominates the orthogonal electric-magnetic 
coherency will be high (assuming the Cagniard theory is 
applicable). On the other hand, if the signal and noise 
are comparable (in one or both series) or if the noise 
predominates in one, then the coherency will be low. 

Figure 16 shows an example of signal and noise of comparable 
magnitude (K = 185-205) while Figure 17 shows a sharp 

peak in the electric at K = 45 with no comparable deviation 
in the magnetic. If the noise predominates in both channels 
then the coherency may be large or small depending on 
whether the noise is common. In our system, coherent noise 
is introduced at the tape recorder. For this reason it 


can be detected not only by the change in ratio but by the 
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fact that it is in-phase except for the shift introduced by 
the analog-to-digital conversion system. For Montreal 

Lake 4-4 H-Ey (Figure 17), noise predominates and is 
coherent for valuesiof K close to 74, 93, 138, 148 and 

Peo To .leo, This has vied to a corresponding wide scatter of 
apparent resistivity values even after decimation on the 


basis of coherency. 


Computational procedures are a second source of 
noise. Although the simple prewhitening filter Zs = X5-Xs_4 
has greatly improved the spectral quality and coherency, 
(Appendix C) some spectra could be improved even further 
by the use of other filters. Figure 18 shows the uncorrected 
power spectra for Hy and E, from Montreal Lake 4-4 using 


data which have been prewhitened by Z, = xX The 


Rage oe! 
Bpecura are far from flat in this region. Even worse 
tne Slopes of the two spectra are quite different and therefore 
> 18 expected that the ratio may deviate considerably from - 
its true value. This example also shows that different 
Drewnitening filters may be required for the E and H 
Sepecora. For notsyirecords.if.has_been, seen that the value 


of coherency using the Daniell estimate is more consistent 


and therefore more useful than the hanning estimate. 


In complex geologic situations the simple linear 
relation between orthogonal components of the electric and 


magnetic field will no longer hold. Under these circumstances 
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the tangential fields will be related by 


Hy = Yi 574 i) ee 
In this case the coherency between orthogonal components of 
the electric and magnetic field may be quite low even in 
the absence of system and computational noise. This situation 
should be suspected if low coherencies are obtained when 


the signal level is high. 


Quon (1963) has shown that for dipoles at iono- 
BonerPic heights, the angle between E and # at the surface 
Of tne earth can deviate considerably from Go", For 


multiple sources this again would lead to low coherencies. 


Prom CHLBY analysis,°it- is’ seen that eliminating 
Gata on the basis of coherency does not discard all the 
data that we would wish and it may discard data from 
complex geologic situations and complicated source systems 


which are meaningful. 


Resistivity Analysis 

im-en attempt to obtain the best. possible Montreal 
Lake apparent resistivity curve for each component, the 
spectral points for which the orthogonal E-H coherency- 
squared is greater than 0.8 have been averaged at each 


frequency. In all cases the high frequency data used was 
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prewhitened using the operator Zs = x Figures 


i ee 
19 and 20 show the resulting apparent resistivity curves. 

Since the curves are quite different in the two directions, 

a layered earth interpretation is not valid. However, to gaina 
semi - quantitative picture, the curves have been fitted 

with two-layer theoretical models. The slope in the 

region 0.2 to 0.4 sec is impossible to fit using a layered 

earth. No attempt has been made to introduce a trend of 

this nature in the theoretical curves as it does not appear 

on the smooth apparent resistivity curves from Montreal Lake 2, 
(Figure o1) but is introduced by the wide scatter of points from 
other curves (Figure 22) which are due to coherent noise 

(Figure 23). It should also be noted that average ae 
apparent resistivity is larger than that of the high energy 


Montreal Lake 2 spectra (Figure 16) for T= .4 to 1.5 sec. 


Averaged phase difference curves have also been 
calculated (Figure 24). A two-layer interpretation of the 
H,-EY curve indicates that the ratio P1/Po lies between 
50 and 100. Using 60 ohm-m as the resistivity of the upper 
layer, the depth to the interface is found to be 17 km. 
Although the interpretation of the H,-E,, curve is very 
uncertain, the ratio p,/p, has been calculated to be 2/1 
with the boundary at 22 km for a 32 ohm-m surface layer. 
The H,-EY phase interpretation is not consistent with the 
apparent resistivity curve. However the H-Ey interpretations 
are in good agreement. Such inconsistencies are expected 


if a layered earth interpretation is applied in an 
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inhomogeneous situation. 


The most striking feature of the resistivity interpretation 
is the apparent geologic anisotropy. Since the magnetic 
declination,.at this station is 1 KZ, the angle between 

the measuring axes and the regional aeromagnetic trend 

is less than 20°. If the aeromagnetic "grain" is 
Poeompanved. bic. Salud lar. charactenistic.<.n wrustalwreaistivity, 
this apparent anisotropy would be explained. Cantwell (1960) 
has. outlined procedures,.to determine. the, major axes..and the 
PesleuLvyioles. along. the. major axes for an, anigeotroplc) earth. 
ie Line case ofan elliptically.polarized wave, this requires 
that the amplitude and phase relationships between all 
horizontal components be known. Since these had only been 
calculated for orthogonal components, the theory using a 
linearly polarized wave was applied even though the records 
show that the field is elliptically polarized. This 

resulted in a meaningless scatter of P. values. The.only 
comforting feature of this experiment was that at high 
frequencies the rotation angle remained close to OF, 


Clearly, too simple a model has been used. 


An attempt has been made to explain the experimental 
data using two-dimensional structures. D'Erceville and Kunetz 
(1962) have considered the problem of two media in contact 
along a vertical fault overlying a horizontal basement that 


is either infinitely conductive, infinitely resistive, or 
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at infinite depth for the electric field perpendicular to 
strike (Figure 25). The only promising model is that of a 
conductive basement with the observing station on the 
conductive side of the fault (P3<P,<Po) and the x-direction 
parallel to strike. Since theoretical data are only available 
at 3 frequencies and for 2 resistivity contrasts (P5/P4 = 9 
and Po/P4 = 100), only a semi-quantitative fit can be obtained. 
Theoretical curves for such a fit are shown in Figure 26. 
In the real situation.the resistivity of the basement will 
De tinive which will decrease the thickness of the upper 
layer and cause a decrease in the phase angle at the longer 
periods as in the experimental data. No theory is 
available for the situation with E parallel to strike 
(H,-E,). However the curves of Neves (1957) for the infinite 
fault indicate that the change in the phase curves will not 
be as largenor the changes in the apparent resistivity curves 
as sharp. This agrees with the observations. A 
major difficulty may be to find a region on the theoretical 
curves where the apparent resistivities differ by almost a 
constant factor over such a wide frequency range. Clearly 
geographic coverage and more theoretical curves are required 
before a final interpretation can be made. It may be found 
that the geology is too complicated to fit with such simple 


models. 
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Observing 
station P 





Figure 25. Finite fault considered by d'tErceville 
and Kunetz (1962). 





ae 
aiitieeuqdi"h yd berslitanges Pie inks 
| (3002) vegan 





mR 
\O 


‘eqjep aye] TeedquoW 4TJ OF peqsnfpe saoeqowered YITM (exTIYS 04 


~~ 


qe[notpuedied |) Jusweseq eATJonpuod YATM 4INeJ TOF SaaAdno [eoTJetosayL “9g sINnSTW 


(09S) potteg (09S) potteg 


OOT OL OOT OT 
Of 


06 os 
seb) 
n 
(4) 
o 
(40) 
0a 

09 6 
O 
n 

OL, 


08 





OT 


OOT 


( W-wWYo ) AQTATISTSer quorzeddy 


— 


my < 


=f 
OGL 


YS 
Cc? al 


wa, 
ial 
~ 


setae] 


mo 

4 5 

< 

) 

"a 

=— | 

- + 

| 

ci 

iD 

ty 

* 

ae = 

= 

= 1 
bow 

oc 

mow 

 } 

at 

ib i 

SF 

@ 

- 

mf 

ma OO) 
- 

“= ade 

; fr 

- 
™ 
* 
pe 
i 
=| 
me 
° 
ede” 
Ve 
tal 
heap 
+ | 


|. 
; | 
; jf 
C yf 
| f 
' 
] 
\ 
\ 


—-y~ 


seiedinaiee ee 


PORERE LE 


> » A reg 
pa Ye a Se 7 4 


foo 


a 
- 


vs 





68. 


In the above interpretation, source effects have 
been neglected. Calculations for a homogeneous half-space 
show that deviations from Cagniard's theory would be unim- 
portant for the conductivities and frequencies under consider- 
ation unless the sources are comparable to dipoles. However 
Tor UiLs Ccompileaved geologic situation and high geomagnetic 
latitude where sources may be particularly complicated, 
GChese etfects may play an important role at longer periods. 
io 1s hoped that by averaging spectra over all records that 
these effects might be "smoothed out." Theoretical 
calculations (Price, 1962; Quon, 1963) show that it is 
Unlikely that the differences at high frequencies could 
result from sources. Further attention will be given to 


these problems in later sections. 


In contrast to the Montreal Lake results, the 
Vulean curves may be fitted by a simple layered model. 
The experimental and theoretical apparent resistivity 
curves, Tor HB are shown in Figure 27. The er 
curves are the same. Sedimentary resistivities were cal- 
culated using logs from nearby wells. Using these values 
for the sediments, the theoretical curves were fitted to 
the observed curves by choosing appropriate basement 
resistivities. The thickness of the 5 ohm-m layer which 
was present at the bottom of the holes has been adjusted to 
give a good fit. This interpretation is basically a three- 


layer conductive-resistive-conductive situation with boundaries 


Bo 3.70 and 35 km. 
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Since only two-layer phase curves are available, a 
detailed interpretation using phase could not be made, 
However, the curves (Figure 28) show the essential features 
of a three-layer conductive-resistive-conductive model. 
Using a two-layer approximation for the upper layers, it is 
found that Po/P 1 ~ 50. For an effective P, = 17.7 ohm-m, 
then p, ~ 885 ohm-m and h, ~ 3.6 km which agrees quite 
well with the apparent resistivity curves. Since the third 
dayer would affect the resistivity ratio, it. is expected 
that Po is somewhat larger. Using phase little can be 
said about 2 except that the contrast would not appear 


bO be nearly as large. 


Roughness of Apparent Resistivity Curves 


Even after decimation on the basis of coherency 
and elimination of those points where coherent noise sources 
are important, deviations in the apparent resistivity as 
large as 100% occur. A second criterion which may be 
used to eliminate points is the ratio of energy in the 
vertical and horizontal magnetic components. For plane 
wave sources and the frequency range of interest, Dosso (1962) 
has shown for a homogeneous half-space that the maximum 
value of | Z/H| ~.01 . In complex geologic situations 
(27H | may become appreciable. For a vertical fault with a 
resistivity contrast across the fault of 1000, Weaver (1963) 


suggests that close to the fault | Z/H | may exceed unity. 
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Quon (1963) shows that for a horizontal magnetic dipole 

over a layered earth Z may become larger than H at long 
periods. In this case the apparent resistivity is much 

less than the true resistivity. Therefore vertical-horizontal 
energy ratios which are larger than normal may either indicate 
a particular polarization of the wave for which the geology 
causes a large He ora source Of small dimensions. - in both 
cases deviations from smooth apparent resistivities can be 


expected. 


In aeSeneral Waywite seems possible to explain 
rough apparent resistivity curves using coherency and 
vertical-horizontal energy ratios. However, in detail, the 
ieoolen 2S more complex. To illustrate this, the apparent 
resistivity and coherency-squared for one set of orthogonal 
components and the vertical-horizontal energy ratio have been 
plotted for two record sections (Figures 29 and 30). In 
addition a smoothed apparent resistivity curve using all 
Montreal Lake records of this component has been plotted 


for comparison. 


The characteristics of the Montreal Lake e-1 
a components (Figure 29) are rather typical of the records 
taken at this location -- the coherency-squared remains 
between 0.8 and 1.0 except for several frequencies where it 
drops to smaller values and the vertical energy is 10 to 


15 percent of the horizontal except for one or two frequencies 
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where it is much larger. A detailed examination of the 
deviations of the apparent resistivity from the smoothed 


mean curve is given in Table 6. It is seen that the maximum 


Table 6. Apparent Resistivity Deviations-Montreal Lake 2-1 H stl 
ee 


Frequency Remarks 
(ops) 
~005 The vertical energy is large and at 
f = .004 the coherency is low. 
.006-.008 Using the criterion that R X R should be 
greater than 0.8, the coherency is marginal 
~O10-.012 The vertical-horizontal energy ratio is large. 
eof =,0L5 Coherency is low. 
Seu Peak in energy ratio. 
.028 Small peak in energy ratio. 
.033-.034 Redes i semarad ned i 
~039 No obvious relationship between the large 
061-.064 deviation, and either the coherency or 


vertical-horizontal energy ratio. Coherent 
noise does not seem to be important. 


a 


deviations are sometimes shifted from the minima of RXR 
and maxima of the vertical-horizontal energy ratio. Also 

in the two cases where the deviation cannot be explained by 
low coherency or high vertical-horizontal energy ratios, the 
coherency is changing quite rapidly. It would seem that 

the windows are not sharp enough to handle the spectral 
changes. However examination of the raw spectra indicates 
that error from this source should be small. It is also 
seen that large deviations may occur in regions where 


the coherency-squared is slightly larger than oo. This 
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confirms the results of Appendix C which indicates that a 


coherency-squared greater than 0.8 is desirable. 


Tape Montreal Lake 5 has very large vertical- 
horizontal energy ratios (Figure 30). With one exception 
the apparent resistivity deviations are highly correlated 
with fluctuations in coherency and the vertical-horizontal 
BPrerey ratio, SAte f = .009 to .O0ll cps” only a small dip 
in coherency accompanies the low apparent resistivity. 


Again this jcoherency-squared is less than. 0.9. 


Baurohner Andlysis of H. and Source Effects at Montreal Lake 


The large vertical-horizontal power ratios of the 
magnetic field at Montreal Lake (Figure 31) compared to 
those in southern Alberta (Figure 32) have led to further 
study of Ho and its relationship to geology and sources. 
Except where otherwise indicated, the following discussion 


is limited to the frequency range 0.001 to 0.070 cps. 


iti .1Llerpreving H at Montreal Lake, the following 
experimental evidence is considered. 
(1) The Hy and A spectra may be quite different in character 


at both Montreal Lake and Vulcan (Figure 33). 


(2) The nH coherency is quite high. Further, it is 


generally much larger than the H,-H, coherency (Figure 34). 


(3) The E, =H, coherency is larger than the HH, coherency 


(Figure 35). 
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(4) The H, spectrum is flatter than the H, spectrum. In 


slope, the ip Spectrum is closest to A (Figure 33), 


(5) Figure 36 shows the power ratio HH in cases where 


the eee coherency-squared is greater than 0.6. For 


Montreal Lake 1 and Montreal Lake 3 the plots follow 


the same general pattern. The Montreal Lake 4-4 


ratios decrease with increasing frequency but are 


somewhat higher. 


(6) Where the H.-H, coherency is low, the Oa ole coherency 


often increases (Figure 37). In this situation the 


H/Hy power Yauio may fluctuate over a rather wide 


range (Figure 38). 


(7) Both ie and H.-H, phases change through 180° in 


the period range T = 2 to 1000 sec (Figure 39). 


(8) For individual events the vertical-horizontal power 


ratio may be much larger than unity (Figure 40). For 


the event marked with an arrow, the ratio is 1.8. For 


the first event in the wave train, Ho is also much 


larger than the horizontal components. 


In a previous section the apparent resistivity 


curves were interpreted in terms of geology rather than 


sources. 


On the basis of the theoretical information available, 


the data may be interpreted as indicating that the field 
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station is on the conductive side of a vertical fault 
striking in the x-direction with a more conductive basement 
at a depth of 10 to 25 km. A large Hoe associated with 

the horizontal gradient of current density component parallel 
to the fault, may be expected. However RokityanskiY (1963) 
suggests, that in the case of the sea-land boundary, a large 
He should only be observed on the resistive side of the 
boundary. This does not agree too well with the calculations 
of Weaver (1963). On the basis of the above observations an 
attempt is made to gain further understanding of the geologic 
eruuacLon.and. ito. determine whether, source effects. also 


contribute to Hoe 


At Vulcan the orthogonal apparent resistivities 
are the same and are quite smooth. In this location, it 
appears that the difference in spectral character is due to 
directional properties of the source. At Montreal Lake the 
geology is more complex. Assuming the source can have 
directional properties similar to those at Vulcan, the 
differences in spectra are likely due to geology as well as 


to directional sources. 


Seeker He near a two-dimensional structure is due 
to the horizontal gradient of current density parallel to 
the structure, and the current density is linearly related 
to Hoe then a high Hie coherence is to be expected. 


Alternatively, if the Hy field were due to a source approximating 
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Beeane Current in the x-direction, and if HY were due toa 
source of somewhat larger dimensions, then a high HH 


coherency compared to the aes coherency is to be expected. 


The fact that the ena a coherency is larger than 
the HoH, coherency shows that H, is more closely related 
Sa .cne currents, than. .to a This indicates that the 
differences arise due to differences in current density 
(geological )rather than to non-plane wave sources. An 
alternative explanation is that the noise level is higher in 
the He system than it is in the EY system. However, at 
higher frequencies, the spectra show that the electric noise 


level is actually the higher. 


since the spectra of Hy and Hi both have the same 
slope at the southern Alberta stations, for purposes of 
discussion it is assumed that they would also have the same 
slope at Montreal Lake were it not for the more complex 
geology or a finite-sized source. Although the calculations 
of Weaver (1963) for E tangential to the fault are not 
in the same frequency range, they indicate that the H/F 
ratio should increase with increasing period. Quon (1963) 
shows that for dipole sources, | Z/H | should increase by 
almost an order of magnitude between T = 100 and 1000 sec. 
For a rather broad north-south current, this change would be 
less. In both cases, an increase in H will be reflected 


by a decrease in thes 
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It has been found that the first four observations 
can be explained on the basis of a two-dimensional structure 
Wien Strike in the x-direction....gim this situation, where 
the Hy -H, coherencies are high, the HL /Hy ratio as a 
function of frequency should be the same on all records. 
Looking only at Montreal Lake 1 and Montreal Lake 3 
(Figure 36), the data might be interpreted-it’ this way. 
However the data from Montreal Lake 4-4 show that the 
relationships are not this simple. The ratios from 
Montreal Lake 5 are even higher although they again show 
a decrease with increasing frequency. These observations 
are probably complicated by the measuring axes being not 
exaculy parallel and perpendicular to strike. Nevertheless, 
iG appears that either inhomogeneities are present for which 
A depends onthe’ directional properties of the source or 


that He changes in time due to a changing source. 


In cases where the ae coherency is low and thus 
little relation exists between them, the HL / Hy ratio should 
be expected to fluctuate considerably. Decreases in 
HH coherency are sometimes accompanied by smaller 
decreases in H,-E,, coherency (Figure 41). This indicates 
that the signal at this frequency may be rather close to 
the noise level. What H, is present may be largely due to 
inhomogeneities which affect the E.. currents. In this situation, 
the H as coherency will be the most sensitive to decreases 


Bg 
in Hy» This would also explain the increase in the H.-H 
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coherency. It also appears that this can be explained by 
an Hy. source of small dimensions but larger than the i 


source. 


No theoretical results are available for the 
phase difference between H, and other components of the 
electromagnetic field close to a two-dimensional structure. 
Quon (1963) shows that for a two-layered earth (p,=5 ohm-m, 
Po = 500 ohm-m, h = 5 km) that (E-H, ) decreases from 
226° at T = 1 sec to 93° at T = 1000 sec for a horizontal 
dipole. However (H-H, ) only changes through 1, None 
of the other dipoles considered provides sufficient phase 
change. It seems that a geological explanation is most 


likely. 


When the vertical-horizontal power ratios for 
several events on a record are much larger than the ratio 
for the entire record, the large ratios for these events 


are most likely due to sources. 


In summary, it seems that all observations involving 
a length of record can be explained on a geologic basis but 
that source effects are required to explain the large 
vertical-horizontal power ratios for individual events. 
Most observations may also be explained on the basis of 
sources but it is unlikely that sources can be used to 
explain the large changes in 9 (E,-H, ) and (H,-H,) or the 


fact that the E_.-H, coherency is larger than the Hy -H, 
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coherency. It may be that over a length of record, geological 
effects predominate; but source effects are important for 
individual events. Clearly, source problems could be better 


Bpoudiedidiny al Simpler’ geologic setting. 


K -Spectral Relationships 
Fp TPE ce Res Senet Onent ps 


Table 2 shows that the ca during the recording 
period at Montreal Lake were of moderate to moderately large 
values. The Montreal Lake records may be conveniently 
divided into 2 groups - those for which ne suG¢wor léss, 
those for which Ke = 44+ , For the records with lower Ky 
the energy falls off quite rapidly up to f = 0.14 cps and 
then decreases very slowly over the frequency range to 
5 cps (Figure 23). The spectra of records with the 
larger Ko decrease much less rapidly to about 
0.75 cps but to a lower energy level (Figure 17). Thus 
it appears that for large — there is a shift of energy to 
longer periods. The one exception to this picture is 
Montreal Lake 4-4 (K,, = 3+) where the spectra has the 
appearance of those with high Bo However for the periods 
on either side of the recording interval K = 4+, It 


Should also be remembered that KY is a world-wide index. 


The three high frequency spectra which were 
calculated for southern Alberta field stations do not show 


such a simple relationship. Vulcan 2-4 (K, = 50) conforms 
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to the above analysis by dropping off quite slowly (Figure 42). 
However the spectra from Brooks 2-2 (K, aH) tar) aot 
very rapidly (Figure 43) while the spectra for Beiseker 1-2 


Be es 4o) are quite different in shape with a broad spectral 


p 
peak in the 3 to 5 sec period range (Figure 44). 


Clearly more data are required to determine whether 
a relationship between Ky and spectra as indicated by the 
8 Montreal Lake records exists. Further, it does not 
appear that there is any direct relationship between Bs 


ana the total energy in this frequency rangex 
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Chapter IV 


SUGGESTIONS FOR ANALYSIS PROCEDURES AND FURTHER WORK 


It has been shown that marked improvements can 
be obtained in spectral quality and apparent resistivity 
values by using a very simple prewhitening filter. It is 
clear that in some cases further improvement could be made 
by using other filters. An ideal way to handle this 
problem would be to have the analysis program calculate the 
Sspectrum,using, unprewhitened. data and fit.a filter to this 
spectrum. . it. should then prewhiten the time series, cal- 
culate its. spectrum, and finally calculate the spectrum of 
the original time series using the power transfer function. 
In cross-spectral analysis the calculated phase shift will 
Beo nave to be adjusted by the difference in phase shifts 
introduced by the filters. With a high speed computer, this 


procedure is feasible. 


The coherency as determined by the Daniell 
estimate has been found to be more meaningful than that of 
hanning but the computational noise level is rather high. 
It is suspected that this is not a property of the Daniell 
window but due to the approximation techniques used in 
CROSP2 . If the noise can be reduced, some prewhitening 
can be avoided. (The speed advantage gained by rounding 
the data before calculating the correlations is as large 


as 15 and so this may not be practical). 
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All data channels should be digitized simultaneously. 
This would minimize the human errors in the conversion 
procedure and allow the computation of phase and coherency 
between all channels with only one conversion per channel. 
The methods of Cantwell (1960) for analyzing anisotropies 
and inhomogenieties could then be applied. It would also 
be preferable to sample all channels at the same instant 
in order that there should be no phase shift between time 
peries. .Kotation of the time series to find the major axes 
would then be practical. Bostick and Smith (1962) found 
that this technique yields better results at the Austin 


station than the procedures suggested by Cantwell (1960). 


The trial and error method of fitting Cagniard curves 
in multilayer cases is time consuming. Moreover one is 
uncertain whether a better fit could be obtained. An 
analysis for automatic curve fitting should be made and 
programmed for the computer. To allow several layers for 
known resistivities of the sediments and a three-layer 


basement, a 5- or 6- layer scheme is required. 


Effort should be made to extend the present theory 
for two-dimensional structures to include the case of 
magnetic polarization and a basement of finite resistivity. 
More type-curves for the present solutions in terms apa 


apparent resistivities are required, particularly as a 
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function of frequency. Additional results from model 


experiments would also be helpful for interpretation. 


The Montreal Lake model can only be regarded as 
a tentative interpretation until geographic coverage is 
obtained. It will be of interest to determine whether the 
directional properties are caused by a single large 
scale discontinuity or by a number of small discontinuities 
miecie GOLreccron of the geologic trend. Experience gained 
Polis area where there is some geologic control should 
be helpful in situations where little geologic or geophysical 


information is available. 


A marked feature of the resistivity interpretations 
eeecnes VarL.arionsin,.crustal.resistivity in. the 3- areas. 
In central Alberta the basement resistivity is about 
1000 ohm-m with no evidence of a decrease with depth 
While.at.Vulcan itis found that the resistivity decreases 
bo about 30 ohm-m at a depth of 35 km. In central Saskat- 
chewan the near surface basement resistivity is in the 
range 30 to 60 ohm-m with a further decrease in the 10 to 25 km 
range. Magnetotelluric techniques are now well enough 
developed that mapping could be carried out to determine 
whether these changes are associated with discontinuities or 


are gradual. 
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In most geologic situations the resolution of 
the magnetotelluric method is insufficient to give any 
detail of the sedimentary resistivities since the product of 
conductivity and thickness of the layers must be large to 
influence the curves. As a tool for investigation of the 
Upper mantley two problems arise. .With induction coils as 
magnetic detectors, locations must be found where the effective 
PeetovL.Viuy OF the sedimenvs,.and crust are nigh in_order 
thac the mantte can be detected with periods of 1000 sec 
Or less. secondly, source effects may become important at 
these long periods. The magnetotelluric method is most 
useful in measuring crustal resistivities. A large scale 
program which should be considered is the mapping of 
crustal resistivities over Wide areas. This would not only 
add to our knowledge of the crust, but also determine 
locations where the mantle can be effectively probed using 


the magnetotelluric method. 
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Appendix A 


Basic Magnetotelluric Theory 


In this section the basic magnetotelluric equation 
is developed following the method of Garland (1960). Fora 
more physical development, the reader is referred to the 


original paper of Cagniard (1953). 


In order to develop the basic equation the following 
assumptions are made: 
oy Micropulsations are plane electromagnetic waves. 
(II) The angle of incidence is zero. 
(III) The earth is a homogeneous half-space. 
DLV) EH is entirely in tne x-direction, i.e., we have only 


Hand H components. 
x iy, 


For an isotropic medium, the time independent 


field equation (MKS units) is 


VE + ‘ate eel rs) eee (1) 


In the earth, for the frequencies under consideration 


oO 
We >> 1 


Hence (1) may be written 


VE + i®oE = 0 (2) 
Under the above assumptions this becomes 
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Yes 
EY = Ae <0 


Since EY mes remain finiteafors Zosiewcg: A) =O 


From (4) and (5) 


s in/4 
(6) 


Therefore E leads H by 45°, 


= 1, (6) becomes 


Changing to practical units and setting w 


2 2 
p= .eT | EL 1/14, | 


Cantwell (1960) shows that this relation holds for 


orthogonal components of # and H in the case of elliptic 


polarization. 


The general propagation constant is 


Ke = UCL + wey 
















i 6 
x 


(4) A Hur - = = 
ei (€) to aotduloa en 
> ve 
eset tt | “Dar 
s(b+ 1 cr = a(t + £) ara} | Lea) 
(2) as a ~ 9A = 2 


x 


wi ae ni 
Oo—= A ..@« & rot svintit oftemet taunt xe eont 


: 7. 
(2) bas () eee 


aN f > 


> (ae) -¥y 


. oe yé H sbssl Z soto . 


‘ 


( a ) 


asmpoesd (6) .[ = 4 artléter bos ation Lsotvoang oF Gamage 
¢ 


y vy 
¥ 
a 


i wi\ 1 81 ms, = 9 


~0o% ebfon nolissie: aids Jsnit gwoda (odeL) Liewinss 
: ney | The) 
obgaiif{e to sesso edd al Ho Bre ai ‘to avaetogmos | -, 


aol LIB LL 
An 


| af Insvenoo aptispegoig [steneg edt id | 
. om : 


i 


ee 4 tuum « 


Pay che ar 


GO =O, therefore k,* = en, 
In the earth 
2 : 2 
Wve << Dol, therefore Ke = Ooi. 


For the frequencies and conductivities under consideration, 
| k,* | S| cae | » hence /sSnell'ts” Law shows that the angle of 
Per,ractiuon is very close’ to-zéro regardless of the angle 
of incidence. Dosso (1962) has calculated the ratio of 
nes amnplilapude=o1eche--vertical-to-horizgontal magnetic field 
for several stratified conductors with geologic parameters. 


The largest ratio found was ~ 14. 


In the more general case of a layered earth, the 
ratio | E/H | can be found by applying the electromagnetic 


poundary conditions. in this .casesthe quantity 


es 2 
eri El /V EN 


pecomes an apparent resistivity (p,) anc bie phase clitererice 
deviates from 45°. Figures A-1 and A-2 show apparent 


resistivity and phase curves for a e-layered earth. 
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Figure A-l. Cagniard two-layer apparent resistivity 
curves. Numbers on the curves show the resistivity 


of the lower medium in ohm-nm. 


The upper medium always 


has a resistivity of one ohm-m and is one km in 
thickness (from Cagniard, 1953). 
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Appendix B 


Spectral Analysis 


The purpose of this section is to present the 
general formulas of spectral analysis, the properties of the 


spectral windows used, and the fundamentals of prewhitening. 


(a) Computational Formulas 


Let Reo tee ee K,. and Tipe #2 ao We 
represent two series of data sampled at time intervals At 
whose means are zero. The spectra of these series are then 
only defined to the Nyquist frequency, fy =1/2At, the 
fuener Irequencies being"aliased into frequencies 

wk Py weer) De. a,.correlation lag index, 

fee, t1,3te).=.4% .+(n-1),Pand) J Ubeva®frequéncy index 
for the (mt+l) equally spaced spectral estimates in the 


ranezener< D,.<¢-f The following quantities are computed. 


N° 
1. Autocorrelations ee ee (vied 
n-L 
1 cna 
AUG) os a X; Xap, ,» A(-L) = A(L) 
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2. Cross-correlations L = -(n-1), (n-1) 
1 n-L 
ee eo) | ea et ie? 
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n-|L| 
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i=l 


3. Autopower density spectra 


n-l 


= 0, (2) 


jean 6 =a ig) 


K(J) = At 


i 


n-l 


Spt) 


Le-(n-1) 


Y(J) At 


where D, (L). is a weighting function - 


(See part (b)) 
4, Cross—power spectrum 


Real part 
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Imaginary part 
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5. Phase lead of the y series over the x series 
Pos OR Pie. BpOye, fp 


W(d ) 


e) = 
(J) Are tan f 7 


with @ taken between O° and 180° if W is 


positive, and between lege and 360° is..W..is negative. 


6. Coefficient of linear coherency 


2 2 
R(g) = V2 J) + W(d 


J)Y( JI 


The autocorrelation function contains the 


emplitude -information- of the spectrum but discards all phase 


Pntormation..-In contrast, the cross-correlation retains the 


phase differences between the two time series. If the xX, 'S 
are considered to be voltages across a pure resistance of 
one ohm, then X(J) Af is an estimate of the power 


dissipated in the frequency interval 


(I-35) ar <r < (3 + 5) af. 
LPothe x, 1S are considered to be voltage and the 


y,'s to be current, then Z-(S) + alee Af is an estimate 


of the power dissipated in the frequency interval. The 


coherency is a measure of how rapidly the cross power phase 


is changing. 
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(bo) The Hanning Window 


In the above spectral estimates the autocorrelation 
functions are modified by a weighting function D, (L) 
called the lag window (Blackman and Tukey, 1958). This 
Meoewne eiltrect of convolving the Fourier transform of 
D, (t) (Q, (f) - the spectral window) with the raw spectra. 


A simple and widely used lag window (whose use is called 


"hanning") is given in continuous form by 


NO) fe 


TT 
Dy (tT) = (1 + cos T) ia |< sue 


= 0 et oe 
m 


where T,, = maximum lag calculated. D,(t) and the 
corresponding Oak iot,.) are shown in Figure B-l. This 
maencdow has the desirable properties of concentrating the 
main lobe of Scat) near f, and side Jobes that fall 
Preto, For digital calculations Dn = ii. OSLTice 
D5(t) =O for |t|> mAt, only m autocorrelations are 


required and the window is computationally inexpensive. 
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Figure B-l1. The lag and spectral window for the hanning 


estimate (after Blackman and Tukey, 1958). 
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(c) The Daniell Window 
The Daniell window (Hannan, 1960; Simpson 
2m ee 1961) exhibits behaviour which can BEaprosach that. or 
the ideal rectangular window. This method uses the complete 
transient (Wiener) autocorrelation of the time function 
X59 tig ee «oh 3 UN 
N-| L | 


1 
Mi eS ee eee e021) ae WR 
fed 


The Daniell estimate is then 





pie sin 2 
fe mm LIT) 5L 
My 2A A(L) Sas ter a) ES 
L=-(N-1) 
ey 
sin— 
where oa is the Daniell lag window. It can be seen 
“ml 


that this estimate has the parameters m and n. In 
Drecvulce 10. ie found._that._che window shape is rather well 
cerinedeby the ratio n/m. Examples of these windows for 
n/m = 10 and 40 are shown in Figures B-2 and B-3. The 
Windows are non-negative; they tend to become rectangular as 


n/m increases; and they are essentially non-oscillatory. 


The basic disadvantage of the Daniell estimate is 
that all lags must be calculated. The window is therefore 


computationally expensive. 
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(d) Prewhitening 


If energy estimates at a frequency ty are to 
be useful, only frequencies close to Ls Should have a 
substantial contribution. If the spectrum rises steeply as 
f leaves Pat an extremely rapid fall-off in the spectral 
window is required. A simple way to avoid this problem 
is to filter the signal before analysis so that the 
spectrum to. be analyzed is fairly constant or slowly varying. 
Since the ideal would be to bring the spectrum close to 
that of white noise, the process is called prewhitening 


(Blackman and Tukey, 1958). 


Maiewcrocedure is simple in principle. Given 


values of x a linear operator is applied to obtain new 


Gb 


Values 2 The relationship between the spectra can 


< 
then be found by calculating the transfer function of 


the operator. 


A simple moving linear combination which is 

@1ven userul is 
Z, = X, + A x, 
The transfer function (Y) is given by the inverse Fourier 


transform of the weighting coefficients. Therefore 


y(f) = 1+ A cos 2rfAt + i A sin enfAt 
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Figure B-4, 


Therefore to obtain improved estimates, the series 
is first prewhitened; the spectrum is calculated; and 
finally, the calculated.spectrum is divided by the power 


transfer function of the prewhitening process. 
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APPENDIX C 


Computational Experiments 


Therpurposesof sthis sectionals torexamine the 
effects on the spectral ratios of simultaneity of conversion, 
trend removal, prewhitening, and different spectral windows. 
Since under the Cagniard theory, the orthogonal components 
of the electric and magnetic field at any frequency are 
Peiearlyerelated, (the usual criterion of quality in» comparing 
signals is coherency. As this relationship is independent 
of time, and the electric-magnetic field ratio is a smooth 
and continuous function of frequency, the smoothness of the 
ratio (or its square) may be used to examine the 
reliability of coherency as an indicator. Although the 
Cagniard theory does not apply at the Montreal Lake station, 
tne interpreted results indicate that the measuring and 
geologic axes nearly coincide and hence the orthogonal 


components are linearly related. 


Effect of Simultaneous Conversion 

The analog-to-digital conversion system used for 
the 1961 data handled only one channel at a time. After 
selecting a record section of interest from the analog 
playback, start and end markers were fixed to the magnetic 
tape. The beginning and end of conversion were triggered 
manually by observing when the markers crossed the playback 


head of the tape recorder. At high playback speeds this is 
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difficult to do accurately and as a consequence we 

invariably ended up with an electric and magnetic time series 
of different lengths having an unknown relative time shift. 
For 3.2 conversions per sec, the normal difference was 

zero to two points. However differences of 18 points have 
occurred. For a tape speeded up by a factor of 16 during 
playback, this difference in length is 90 sec. Manual 
alignment of digital data is slow and unreliable. This leads 
PoeeLoss of contidence in computational results. ‘This 
section considers the effect of these shifts on spectral 
results by comparing independent conversions of given 

analog records. Bogert (1962) compares manual and electronic 
conversions of the same analog records. In a sense, this 


Peowon 18 a-continuation of his work. 


Various parameters can be used to compare two 
records including their variance, auto- and cross-spectra, 
coherency and phase. The simplest of these is variance. 
Prcnough if contains no indication of the way the energy 
is distributed in either frequency or time, it will show 
gross differences between independent conversions. Auto- 
spectra are more sensitive, but a pair of records can have 
identical spectra yet be incoherent. Hither point-by- 
point equality in time, or identical autospectra and unity 
coherency serve as sufficient criteria for identity of a 


pecord pair. 
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The first experiment consisted of analyzing 
slightly shifted segments of one digital record. The 
eross-—spectra of points 1 to 1000 with points 11 to 1010, 
el to 1020, and 31 to 1030 were calculated. The results 
will of course depend on the differences between the 
non-overlapping segments. Figure C-1 shows the spectrum 
of points 1 to 1000; Figure C-2 shows the effect of shift on 
spectral ratio; and Figure C-3 shows the coherencies. 

The percentage differences in variance were 3.5, 6, and 10 
Percent. =Ll- 1S -apparent that small shifts in real 

Breda ts-canaleéad to substantial differences in the ratios. 
PipesuniaLlcaves that a particular spectral bump can often be 
associated with a discrete event on the analog record. This 


Observation is borne out by many other examples. 


Both coherency-squared and the percentage differences 
deteriorate with increasing record shift, as expected. 
There 18 a correlation of large positive percentage difference 
with low coherency and spectrum troughs in each case. 
However the outstanding percentage differences (ignoring 
sign) occur between frequencies of 10 and 20, where 
coherencies are only slightly lower than elsewhere. This 
example points out the typical behaviour of coherency as 
an indicator of ratio smoothness. Large coherency is 
necessary, but in itself is not sufficient to ensure a 


Smooth spectral ratio. 
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In the section on spectral programs is an example 
of the spectra resulting from a two-channel conversion of a 
Single analog signal. Any differences must be due to the 
Pereriuedttrerence in sampling’ instant, and to minor 
Brrrerences in the response and noise of the playback filters. 
Although the vast majority of values of coherency-squared 
Bapccoe. 3, co Lew speccral ratios greater than 1.5:1 do 
emerge. As is to be expected, the major differences occur 
for low spectral amplitudes. Other examples of two-channel 
conversion of a single signal have led to the same conclusion. 
The process is very effective in improving coherency but not 


furvesas (Successful in ratio smoothing. 


One example of differences resulting from 
independent conversions is illustrated by Figure C-4. This 
Shows a pair of spectra derived from separate conversions of, 
ostensibly, the same segment of analog Se eeicne The 
difference in variance is 6 percent. Table C-1 summarizes 
variance differences found in independent conversions of 


e7 different analog records. The number of points were 


between 4000 and 5000 in all cases. 


eters Snr seers SS ae 


i It is unlikely that any appreciable part of the difference 
can be attributed to conversion round off. 
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Table C-l 


Ditrerence «in Variance Number of Cases 

b 

O-1 Le 

els fr. 2 6 

Bele = 3 4 

6 1 

24 2 

38 1 

87 if 


It should be emphasized that near-equality of variance does 


not guarantee near-identity of spectra. 


The problem of ensuring conversion of corresponding 
sections with a single channel analog-to-digital conversion 
System can probably best be solved by post-recording pulses 
Sieaecirect record channel of the tape. These pulses can 


tnen be used to trigger the converter. 


Trend Removal 

With the data acquisition and processing equipment 
used, and the spectra of the signals recorded, our raw 
digital data must be expected to contain finite power at 
zero frequency as well as finite power in the lowest energy 


band. Both of the spectral programs used for the 1962 data, 
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CROSP2 and CS TUKS, remove the mean but neither removes 

the slow trends. The program written here for the IBM 1620 
and used in processing the 1961 data removes a linear trend 
Desed on the means of the first and last thirds of the 
record. Almost invariably the lowest frequency band contains 
Seraree percentage of the total record energy. If it were 
pure signal the spectra would probably be useable although 
the time series would first require prewhitening to prevent 
energy from leaking through the side lobes. More often 

than not, however, a large contribution is due to the baseline 
miaac Of the amplifiers used in the field due to ambient 
temperatures. The importance and techniques of trend 


removal are discussed by Blackman and Tukey (1958). 


The experiments were of two kinds. One consisted 
of computing a number of spectra of real data, both with 
and without linear trends removed to determine the change 
in ratio and coherency. The other was the computation of 
the spectra of segments of sinusoids, for which the power 


Spectra are known analytically. 


GratsinskiY (1962) shows that the spectrum of a 
fraction of a sinusoid (starting from a zero-crossing) 
undergoes a number of systematic changes as this fraction 
increases from 0.5 to 2 or more. For 0.7 of a cycle or less 
the spectral peak is shifted toward zero frequency and the 


energy is spread fairly uniformly over the spectrum up to 
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twice the frequency of a pure sinusoid. Some of his 

examples have been reprogrammed and are shown in Figure C-5, 
The results using the IBM 1620 program with mean and linear 
trend removal are shown in the same figure for two records, 

one consisting of the first 198 points of a 600 point cycle, 
the second consisting of the first 300 points of the same 
Beers i Dov, casés m= 100. All the data were 

normalized to have the proper variance. The results illustrate 
that mean and trend removal significantly decreases this 

long period energy and add energy to the higher frequency 


regions of the spectrum. 


Inasmuch as the drifts actually observed are often 
nearly linear, it should be hoped that the subtraction of 
idinear trend would eliminate a majority of the drift and 
Emprove the spectral estimates. The effect will of course 
eeeoecdepend on-the distribution of true energy in the lowest 
frequency band and the overall spectral shape. In Figure C-6 
are shown two spectra calculated from the same digital record, 
with and without trend removal. In both cases the mean has 
been subtracted from the record. Figure C-7 shows the 
coherencies in the two cases, and Figure C-8 the resultant 
Bpecctral ratios. Included in this last figure is the ratio 
computed after prewhitening the two series, as an indication 
Of what this ratio "should" be. It is evident that the 
trend removal has decreased slightly the power in the first 


two bands but has added power to all others beyond the 
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EFFECT OF TREND REMOVAL 
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third. It is also evident that the amounts added to the 
two spectra (electric and magnetic) were not the same, 
since the ratios are consistently larger in the trend 


corrected E-H pair. 


In this example, removing the trend has reduced 
the coherency and the accuracy of the ratio in almost every 
frequency band including the lowest few. This example may 
be anomalous, as in other cases, trend removal has led to 
higher coherency and smoother ratios. However it does seem 
that trend removal must be used with considerable caution. 


fouseprocedure will be discussed further in the next section. 


Prewhitening 


If the spectrum changes rapidly then prewhitening 
of the data before computing the spectrum is recommended. 
bince the raw high frequency spectra fall off approximately 
as ‘yaa a filter which strongly rejects the low frequencies 
is required. Because of the difficulties of using a computer 
at large distances, the only prewhitening filter used was 


Z The power transfer function of this filter 


meee ret. 
is shown in Figure B-4 (A = -1). This procedure invariably 
gave much improved ratios and coherencies. The low frequency 
(LF) data were not treated, since, with two exceptions, 


their spectra were already relatively flat and their 


coherencies high. (The low frequency data were subjected 
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to trend removal. In view of the results of the previous 
Section, it appears likely that a form of prewhitening 
similar to that used for the HF data, but less drastic, 
would have given better results than the removal of linear 


trends. ) 


Figures C-9 and C-10 show HF spectra and 
coherencies without and with prewhitening, respectively. 
Pigure C-ll shows the resultant ratios in the two cases. 
These ratios are computed from the raw spectra and are 
tilted by the difference in response curves of the electric 
ana Magnetic detection systems. Other examples of the 
effectiveness of prewhitening in smoothing spectral ratios 
can be seen in Figure C-8 and by comparing Figures C-19 


Several sets of HF data were analyzed both with 
and without trend removal prior to prewhitening. It was 
found that the maximum spectral differences were of the 
order of 1-2 %. Thus, if one is going to use prewhitening 
filters which incorporate drastic long period rejection, 
trend removal is normally not required. This is, of course, 


expected. 


Spectral Programs 


Our early analyses included two types of data which 


are very demanding of spectral computations: sinusoidal 
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calibrations and untreated HF data. The digital data were 
obtained by a single channel conversion and the spectra 
calculated using CS TUKS which computes a hanning estimate. 
One source of difficulty arises from the negative side lobes 
of this window. These resulted in negative spectral 
estimates over as much as 90% of the spectrum when applied 
Go the 1/f* HF spectral functions. We also attempted to 
calculate the frequency response of the detection systems 

by comparing the computed magnitudes of spectral peaks of 
each component with those computed for an input calibration 
signal. The results were very erratic. This was probably 
meg an Pare, toOuslipntly dififerent conversion rates 

(< .01%) and to not quite identical record segments. Another 
source of error was the contamination of small peaks by the 


side lobes of large peaks. 


About this time the program CROSP2 (Simpson et al, 

1961; 1962; Galbraith, 1963) became available. Since it 
computes a Daniell estimate and is very fast, it seemed 
ideal for our purposes. Micropulsation "pearl" data, 
previously analyzed using the hanning window, were kindly 
run for us by Prof. T. R. Madden using a Daniell estimate. 

A comparison showed that the Daniell window resolved the 
spectral fine structure of the pearls much more clearly on 


the sides of steep slopes, as shown by Figure C-12. Ratios 
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computed from these Daniell estimates in the same frequency 
range were also more consistent with the known electrical 
properties of the area. Coherency behaved much more reasonably 
in that values were always between zero and one. It is also 
Pound thal they do not change drastically with m. As a 
consequence, for a majority of the 1962 analyses CROSP2 

Mas used. Atlew data sets have since been found for which 

CS TUKS gives better results than CROSP2. (Other spectral 
programs are available, notably BOMM (Bullard ét al, 1962) 


but CROSP2 would appear to be desirable for its speed.) 


The two windows were compared by using the programs 
to analyze several kinds of signals: computer generated 
sinusoids’ cfitvarious amplitudes and frequencies, segments 
Gnicalitbration isienalwdigithzed from *tapes recorded in the 
mee. eandiinatural datas relt should “be spointed outethat 
CROSP2 rounds the records after removing the mean, but 
before computing correlations. The errors introduced will 
not exceed 1 part in 800 of the range of data values. 

GS TUKS uses the full data accuracy. This emphasizes that 
the results of computational experiments with CS TUKS and 
CROSP2 may reflect not only the different window properties 
but also the particular computational procedures involved 

in the program. Accuracy of the data itself is limited by 
the 3 digit plus sign fixed point form (1 part in 2000 of the 


maximum range). Thus at frequencies where the Signal energy 



























esend motlt besvamos 
jelenos sTom dele ‘stew onNBT 
rola .BS18 edt to sstiteqorwg 
Wisd 1 wis 9tew esuls BAS Os 
‘3 gard boavot 

oONsUpSsanes 

gE BISD WS O98 @5wW 

d esvian GHUT Bo 
lieve o%5 ene tRacrees 


‘Scr r ff 


a: Lvuow %S Und 3o@ 


shri] fsrevyse envisns oF 


iiiqmns evoinsv 30 sbiosun 


abmauo 29 208) 
ro» satsvanoo o10tst 
nf wmaq I bessxs 

Sa SJB iiud sas ‘egegu CAUT 
isnoltsecsuqmoo to evivuesY 
‘3 ylao Jon tosliter vem Sia 
quod tsivotiteq sat opie Ba 
iO (YOSILOOA .heTEbia ene 
q bexl! cate auig ghete & 


pst) te auclT .(eznst mumie 
uj 


7 


=. 





is sufficiently low, the CROSP2 round off becomes evident in 
the spectra (Figure C-9) which are overly smooth, coherencies 


which are low, and ratios which are meaningless. 


The two series 


34 (kK)o= 200° sin’ (2kir/?) 
and 
So (ky) =*900"sin €2kr/T) 
feerece i= J 2)... . . 4000. (approximately) 
Va Der) Od 


were analyzed by each program to see whether the correct 
ratio was indicated on the peaks, and how this ratio changed 
with distance from the peaks. In both programs ratios 
computed at the peaks were correct to within one percent. 

an cre CS°TUKS output this same spectral ratio was 

Maintained to within a few percent over the entire frequency 
Panee Of the analysis. Ratios of the two CROSPe spectra 
Becrease moroionically to the order’of unity at higher 
frequencies. The CS TUKS spectra alternated between positive 
and negative values beginning two points away from the peak 
whereas the CROSP2 spectra were all positive with the exception 


of a few points well down into noise. 
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Figure C-13 shows the spectra of S, (k) fora 


2 
series of frequency f, = 061 fy - As f increases beyond 
a both spectra fall off 5 or 6 orders of magnitude with 

CS TUKS decreasing more rapidly with alternate values being 
negative. For f less than ee ChUsracontly, Pail saories 
orders of magnitude to the program noise level while 

CS TUKS is down in excess of 5 orders of magnitude. The 
precise behaviour of the spectra depends on the place in 
the cycle where the series are terminated, and in the case 
fee noote,2 On. Lnputedata format. The changes in. relative 


energy dependent on these factors is negligible. This is 


expected since the series contained over 150 cycles. 


These results indicate that the erratic 
Saltpration results using CS TUKS may not have been the 
Peuilt Ofatnesprosram. [i appears likely that the errors 
may have been caused by comparing slightly different 
BeevLone oi rcape, Weinge slightiy different conversion rates, 
and if not enough cycles were used, different starting and 
PeempactOn points of theicycles. If errors are 
introduced by the side lobes, the problem could be resolved 
by analyzing one frequency at a time. Further experiments 


using multi-channel conversion are recommended. 
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As applied to CROSP2, this problem has been treated 
by analyzing calibration signals converted by the 2-channel 
analog-to-digital conversion system. In the first example, a 
square wave of frequency .0961 cps was converted 13.6 
times per second on two channels whose sampling instant differed 
by 12.4 ms. Auto~ and cross-spectra were computed for the 
first 1000 points with m = 250. One of the spectra, coherency, 
and percentage difference between the two spectra are plotted 
in Figure C-14. The average difference of 2.3% is due to 
insertion loss differences between the two playback filter 
Systems. The short dashed lines over the spectral peaks 
are the heights (relative to the first peak) predicted for 
Saeotifsnloe number of input cycles. With increasing 
jreouency he peaks computed for the odd harmonics diverge 
Prom uneir predicted values but the amount of energy in 
the even harmonics is seen to increase. At frequencies 
greater than 4 cps the playback filter reduces the signal 


at a nominal rate of 24 db per octave. 


There are a number of interesting features in 
this figure. Ratios at odd harmonic peaks fall within 
2-3% of the mean well beyond the Hoes? harmonic. There is an 
increase of the (smoothed average) ratios with increasing 
frequency, to about +10% near 5 cps, beyond which they again 
diminish. This suggests response differences in the two 
playback filters. Despite this, 50 percent of the power 


ratios lie within 6 percent of the mean. 
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haocner point of interest is; that, if ratios 

are discarded on the basis of coherency the scatter becomes 
much less. There appears to be a limit (at R = .95) above 
which this is no longer very effective. More effective is 
the process of also discarding ratios adjacent to points 
where R< .95. Very few ratios near the mean, and all 

but one of the extreme values are eliminated in this manner. 
mince the twovseries used are derived from a single input 
Sven *convercea at two slignhtiy different instants, the 
results are also of Inveres. “TO Our eariver discussion of 


conversion tecnniques. 


In the analysis of real calibration signals, 
spectral ratios versus frequency were computed. Their 
ratios agreed with predicted values to within two percent 
Wise 1S tne accuracy of the predictions themselves. One 
Siimanuaeeuo, OlLeital processing of calibration data is 
thao it serves ae a check of the entire record and playback 
svevens MH "Second advantage lies in the Tact that, in the 
magnetic systems there is always some natural signal 
superposed on the calibration, Digital processing permits 


CGallbrating even during periods of large natural signal. 


ratte 


An example is next shown in which CS TUXS indicates 
larger coherencies than does CROSP2, although the ratios 


are of comparable smoothness. One pair of digital E-H 
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records had mean and trend removed. They were then analyzed 
by CS TUKS and CROSPe. The records were then prewhitened 
and reanalyzed. The spectral pairs and coherencies-squared 
are shown 1n Figures C=-15,°~C-16; C-17, and°Cc-18. The E/H 
ratios are shown in Figures C-19 (not prewhitened) and C-20 


(prewhitened). 


Finally, Figure, C-2l ilillustrates the difficultiés 
of using CS TUKS coherencies as a measure of spectral quality. 
Although the spectra appear to be unrelated, the coherency 
Waties trom near @ero to much) greater’ than one (Ro Pe 1B 
piloeced-at Ao 1). The Daniell window also encounters 
difficulty (Figure C-9). However the CROSP2 coherencies 


become very small indicating the unreliability of spectral 


values. 


Conclusions 

With several qualifications, it is concluded that 
digital power spectrum analysis can be used as a qualitative 
tool in the processing of magnetotelluric data. By operating 
on series whose spectra are relatively flat, the accuracy 
of estimates is significantly enhanced and the range of 
frequencies over which reliable estimates are obtained is 
increased. Prewhitening was the most effective process 
found for smoothing ratios and increasing coherency. 


Removal of a linear trend does not always improve the results. 
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A 50. 


Diagram 44 


Illustration of HF spectra without prewhitening (CROSP2). 


Figure C-17. 
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If possible, undesirable energy should be removed from the 
lowest frequency band by a prewhitening filter rather than 


by subtraction of the trend. 


In the case of spectra which contain one or more 
spikes, Daniell windows were found to be more satisfactory 
than were the hanning, in that the Daniell estimates are 
all positive, all coherencies are in the range zero to one 
(hanning estimates are often greater than one) and the 
spectral ratios are as predicted to the order of one or 


two percent. 


It is necessary to synchronize digital records if 
bney are to be compared. This is obvious for comparison 
in the time domain and our results show that they must be 
erosely aligned if their spectral ratios and coherencies 


are to be meaningful. 


Coherency is only a rough measure of the quality 
of data and spectra, as judged by spectral ratios. Even 
very large values of coherency (> .99) do not ensure smooth 


Tavlos. 


For a given digital record, the spectra and 
consequently coherency depend on the computational techniques 
used in the particular program. Hence, values of coherency 
cannot be directly compared. Some experience is required with 
each program to enable one to judge the significance of its 


coherencies in a particular situation. 
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Appendix D 


PLOTS OF SPECTRAL DATA 


Notation and Conventions 
As an aid in using the graphs, the abbreviated 


notation used on the plots is explained below. 


Notation Meaning 


DELTF - PTrequency increment 4t¢ wnoieh 
parameters are calculated 


FMAX - maximum frequency calculated or 
plotted on graph 


K - frequency index (true frequency at 
K is K times DELTF) 


LAMBDA - 1 - indicates data have been prewhitened 
using Z, = X, - X4_4 

RHO - apparent resistivity 

SS aaa! - coherency squared 

WEIGHTED COHERENCY _ Fy(R, XR.) + PCR, X RL) 


P.. + Ey 


where Pie is are He and H spectral 
densities; 


: 
Z 


4 


R_ ¢ -H| and H_-H 
age hg are HL ' + 


coherencies-squared 


XX - indicates exponentiation 


On some graphs the term "COHERENCY" is used while 
on others "COHERENCY XX 2" appears. In both cases these 


should be read as coherency-squared. 
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Diagram 25 will be used to explain the numbers 
Ba oie recenas! deontiti cationwline,. .".7 "indicates ‘the 
analog-to-digital conversion run while "1 TO 1087" indicates 
the points from this conversion that were used in the 


analysis. 


oe peptendeeta-spéectral plots have been corrected 
por anplitier settings. filtering during conversion, and 
prewhitening (if used). The electric spectra have also been 
eorrected for the response of the detection system while for 
Dne magneticsrit has been plotted on the same graph. 
PeEnoury che graphs indicate’ an absolute level for the 
spectra, there is an insertion loss of 41.5 percent in the 
analog-to-digital conversion system which has not been 
taken into account. Values of coherency-squared which 
are greater than one are plotted as equal to one. On 
graphs which are plotted as a function of K, the calculated 
values are plotted at integer values of K, with points 
between plotted using linear interpolation (except for 
Bpparent resistivity curves for which there is no 
interpolation). On graphs, which are plotted as a function 
of period, there is no interpolation. In this case, those 
points for which T is the same to two significant figures 
are averaged before plotting. All plots of field ratios 


are in terms of power. 
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Arrangement of Plots 


This appendix contains the machine plots used in 
Piperprevacion of the spectral data. For completeness of 
Paci SET, Chose used in the text as examples are repeated 


in this section and are designated by the same diagram number. 


Mie OLOLe are arranged 1n two sevs: 
(a) Plots of apparent resistivity and phase versus period. 
The averaged Montreal Lake plots are followed by 
the plots from individual Montreal Lake tapes and the 


southern Alberta data in alphabetical order. 


(ob) Plots of other parameters. 


The plots from each tape are grouped together. 
The Montreal Lake data are first with the tape number used 
to determine sequence and the southern Alberta data following 
in alphabetical order. The data from each tape are 


arranged as follows: 


Low Frequency (CROSP2 used for all spectral calculations) 


(1) Hy and EL spectra, coherency- squared, Hy calibration 


(2) Hy and BY spectra, coherency-squared, Hh calibration 
(3) He spectrum, vertical-horizontal power ratio, weighted 
coherency-squared, H, calibration 


(4) H.-H. and H.-H, coherencies-squared 
> aie 7 aaa? 
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(5) H,/H, and HL / Hy power ratios 
(6) apparent resistivity versus frequency (southern Alberta 
data only) 


(7) any additional plots 


High Frequency 


(8) Hy and BS spectra, coherency squared, and Hy calibration 
(CROSP2 , prewhitened) 

(9) A and EY spectra, coherency squared, and A calibration 
(CROSP2 , prewhitened) 


(10) any additional plots 
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